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ABSTRACT 


A  brief  summary  of  Che  results  obtained  in  research  sponsored  by  the 
Naval  Research  Laboratory  under  Contract  N00014-80-C-0802  is  presented. 

The  research  covered  several  problems  in  the  area  of  spread-spectrum  random- 
access  communications  for  fading  channels.  The  results  are  applicable  to 
the  Navy's  intra-task-force  communications  network. 
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SUMMARY  OF  RESEARCH  RESULTS 


Our  work  in  spread-spectrum  communications  that  was  supported  by  this 
contract  focused  on  the  performance  of  slow-frequency-hopped  (SFH)  spread- 
spectrum  communications.  A  variety  of  channel  models  were  considered  to 
reflect  varying  degrees  of  amplitude  fading  and  selectivity.  Of  primary 
interest  is  the  class  of  slow,  nonselective,  Rician  fading  channels  in 
which  there  are  two  components  of  the  received  signal:  a  non-faded  or 
direct-path  component  and  a  faded  or  scatter  component.  The  scatter  com¬ 
ponent  is  assumed  to  undergo  slow,  nonselective,  Rayleigh  fading  in  this 
model,  so  that  the  Rayleigh  fading  channel  is  obtained  as  a  special  case 
(no  direct-path  component).  In  addition  to  the  received  signal,  additive 
white  Gaussian  noise  is  present  at  the  front  end  of  the  receiver,  so  the 
non-faded  additive  white  Gaussian  noise  channel  is  also  obtained  as  a 
special  case  of  the  general  model  (no  scatter  component).  The  effects  of 
selective  fading  were  also  considered,  especially  frequency-selectivity 
which  produces  intersymbol  Interference.  Details  of  the  various  channel 
models  are  described  in  Appendix  A. 

Several  bounds  and  approximations  for  the  bit  error  probability  in  a 
SFH  spread-spectrum  multiple-access  system  are  presented  in  Appendix  A  and 
in  [l]-[3]  .  Both  FSK  and  DPSK  data  modulation  and  selective  and  nonselec¬ 
tive  fading  channels  are  considered.  These  results  are  very  general  in 
nature  and  can  be  adapted  to  a  wide  range  of  systems  and  channel  models. 

A  specific  problem  that  can  arise  in  a  system  like  the  intra-task- 
force  (1TF)  communications  network  is  due  to  the  possibility  of  a  specular 
multipath  signal  with  a  relative  delay  greater  than  the  dwell  time  of  the 
frequency  hopper.  This  problem,  although  not  addressed  specifically  in 
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Appendix  A,  can  be  analyzed  by  the  results  developed  in  our  research.  The 
results  of  Appendix  A  are  applicable  to  a  system  in  which  there  are  K  simul¬ 
taneous  SFH  signals.  If  there  are  K'  simultaneous  transmitters  and  if  each 
signal  produces  one  nonselective  Rician  faded  component  plus  one  specular 
multipath  component,  then  there  are  K*2K'  interfering  signals.  If  the 
relative  delay  of  the  specular  multipath  component  exceeds  the  dwell  time 
then  the  bounds  and  approximations  given  in  Appendix  A  apply  with  K  * 2K* . 

In  such  a  system  there  are  K-l  *  2K’-1  interfering  signals  for  each  receiver. 

Application  of  the  results  of  Appendix  A  to  such  a  system  is  illustrated 
by  the  data  of  Table  1.  The  only  reason  for  presenting  numerical  results  for 
this  special  situation  is  that  perhaps  larger  values  of  K  are  of  interest 
than  for  a  system  without  the  specular  multipath  components.  The  results  of 
Table  1  are  for  the  same  model  as  described  in  Appendix  A,  and  the  notation 
is  exactly  as  used  in  Table  2  of  Appendix  A.  The  approximation  PG  given  in 
Table  1  is  a  new  result  (obtained  after  [4]  was  submitted  for  publication), 
which  we  believe  to  be  slightly  more  accurate  than  the  approximation 
described  in  Appendix  A.  Both  and  are  for  channels  with  fading 

which  is  slow  relative  to  the  hopping  rate  (case  (i)  described  on  page  10 
of  Appendix  A) . 

The  methods  and  results  developed  in  Appendix  A  can  also  be  applied  to 
determine  the  probability  of  error  in  a  coded  SFH  spread-spectrum  system. 

For  a  fully  interleaved  system  these  results  can  be  applied  directly.  This 
is  because  the  interleaving  breaks  up  the  error  bursts  due  to  the  fading 
and  multiple  access  interference,  in  which  case  the  bit  error  probability 
is  the  performance  measure  of  interest.  Thus,  for  interleaved  systems  the 
performance  of  various  codes  can  be  determined  from  the  results  given  in 
Appendix  A  and  published  data  on  the  performance  of  the  codes  for  the  binary 


memory less  channel. 
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Table  1.  Bit  error  probability  for  nonselective  Rayleigh 
fading  and  specular  multipath. 

a)  K  -  10  (K*  -  5),  q  -  100,  and  Nb  =  10 


S/N0  (dB) 

?L 

PG 

PA 

pu 

6 

0.161 

0.175 

0.182 

0.199 

8 

0.118 

0.131 

0.137 

0.156 

10 

0.085 

0.096 

0.102 

0.123 

12 

0.060 

0.071 

0.076 

0.098 

15 

0.036 

0.046 

0.051 

0.074 

20 

0.018 

0.028 

0.032 

0.056 

OO 

0.009 

0.019 

0.023 

0.047 

K  -  20  (K' 

*  10)  ,  q  ■ 

100,  and  * 

10 

&/N0  (dB) 

?L 

PG 

PA 

pu 

6 

0.153 

0.187 

0.193 

0.230 

8 

0.115 

0.145 

0.151 

0.192 

10 

0.085 

0.113 

0.118 

0.162 

12 

0.063 

0.089 

0.093 

0.140 

15 

0.042 

0.066 

0.069 

0.119 

20 

0.025 

0.049 

0.052' 

0.103 

OO 

0.018 

0.041 

0.043 

0.095 

K  *  50  (K' 

-  25),  q  =■ 

250,  and  - 

10 

2/Nq  (dB) 

?L 

PG 

PA 

pu 

6 

0.157 

0.188 

0.191 

0.227 

8 

0.119 

0.147 

0.149 

0.189 

10 

0.088 

0.114 

0.116 

0.158 

12 

0.066 

0.090 

0.091 

0.136 

15 

0.044 

0.067 

0.068 

0.114 

20 

0.028 

0.050 

0.050 

0.098 

OO 

0.020 

0.043 

0.041 

0.090 

Table  2.  Bit  error  probabilities  for  uncoded  and  coded  SFH 
systems  with  nonselective  Rayleigh  fading. 

(K  -  15,  q  -  1000,  and  1^-40) 

■'«!«•>  pb 

4.0  X  10~?  1.1X10'P 

2.7x10  ,  1.1X10  * 

1.9x10,  1.7X10,- 

1.3  X  10~^  4 .0  X  10~1J 


14 

16 

18 

20 
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For  a  system  without  full  interleaving  the  bit  errors  are  not  indepen¬ 
dent  and  thus  the  bit  error  probability  does  not  completely  describe  the 
channel  performance.  However,  we  have  evaluated  the  performance  of  certain 
Reed-Solomon  codes  with  partial  interleaving,  and  typical  results  are  given 
in  Table  2.  The  approximation  PA  to  the  probability  of  error  for  an  uncoded 
system  is  compared  with  the  bit  error  probability  for  a  system  which  uses 
a  (255,127)  Reed-Solomon  code  with  partial  interleaving.  Notice  that  for 
values  of  ^/Nq  greater  than  18  dB,  the  coded  system  gives  several  orders 
of  magnitude  improvement  in  the  bit  error  rate.  Further  work  on  the  perfor¬ 
mance  of  coded  SFH  spread-spectrum  systems  is  in  progress  (primarily  under 
other  sponsorship). 

The  data  in  Table  2  gives  a  comparison  between  the  performance  of  uncoded 
and  Reed-Solomon  coded  systems.  Another  interesting  comparison  is  the  perfor¬ 
mance  of  a  Reed-Solomon  coded  SFH  system  for  two  different  sets  of  assumptions 
on  the  frequency  hopping  and  interleaving:  (i)  no  frequency  hopping  and  no 
interleaving  vs.  (ii)  frequency  hopping  with  interleaving  of  the  code  symbols. 
The  channel  model  that  we  consider  for  this  comparison  is  the  very  slow, 
nonselective  Rayleigh  fading  channel.  In  case  (i)  we  assume  that  the  instan¬ 
taneous  power  in  the  received  signal  is  constant  for  the  duration  of  the  code 
word,  but  in  case  (ii)  the  instantaneous  power  is  constant  for  the  duration 
of  a  code  symbol  but  (because  of  hopping  and  interleaving)  the  power  levels 
for  different  symbols  in  the  same  code  word  are  independent.  In  Table  3 
numerical  values  for  the  block  error  probability  are  presented  for  the  (31,15) 
and  (255,127)  Reed-Solomon  codes.  The  probabilities  and  are  the 

block  error  probabilities  for  cases  (i)  and  (ii),  respectively.  The  data 
is  for  a  system  with  only  one  transmitter  (K  ■  1)  in  order  to  isolate  the 
effects  of  fading  from  the  effects  of  multiple-access  interference. 


Table  3.  Block  error  probabilities  for  a  coded  SFH 
system  with  nonselective  Rayleigh  fading. 


(31,15)  Reed  Solomon  code 

S/Nfl  (dB) 

p(i) 

E 

p(ii) 

E 

-2 

-5 

20 

8. 29  X  10 

8.44  X  10 

22 

5 . 31 X  10  2 

2.50  X  10~o 

24 

3.39 x  10  j 

5.92x10"° 

26 

30 

2.15  X  10“~ 

8.63  X  10“ 

1.21 X  10  * 

3. 97 x  10 

(255,127) 

Reed  Solomon  code 

S/N0  (dB) 

p(il) 

E 

15 

2.85 X  10“| 

8.86 X  10“^ 

16 

2.34  X  lO-!* 

1.50X  10“^ 

17 

1.90 x  io“: 

3.52X  10  l 

18 

1.55 x  io“: 

1.53 X  10, 

19 

20 

1.25 x  io"r 

1.01 x  10 

1.65 X  10"  3 
5.62 X  10"iB 
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A  significant  area  of  progress  in  the  random-access  area  under  this 
contract  has  been  the  design  and  analysis  of  retransmission  control  poli¬ 
cies  for  a  random-access  broadcast  channel  [51  —  [7],  [9]  .  The  policies  can 
be  implemented  in  a  distributed  fashion.  Analysis  of  delay  and  throughput 
is  provided  in  these  papers  using  the  concept  of  local  Poisson  approxima¬ 
tion  which  is  introduced  in  these  papers. 

Versions  of  the  recursive  retransmission  control  policies  which  are 
relatively  insensitive  to  the  traffic  statistics,  and  modifications  which 
reduce  feedback  information  requirements  are  also  reported  in  [5] . 

It  is  proven  in  [71  that  the  retransmission  policies  in  [5]  provide 
stable  throughput  at  rates  of  up  to  e  ^  packets  per  slot.  Moreover,  a 
general  methodology  for  proving  such  stability  results  is  provided  in  [7] 
and  the  methods  are  also  applied  in  [7]  to  prove  a  strong  stability  property 
of  G/G/l  queues  which  is  of  general  interest  for  queueing  network  studies. 

Even  though  the  papers  [5]-[7],  [9]  do  not  deal  explicitly  with  a 
spread-spectrum  system,  they  were  developed  for  spread-spectrum  applica¬ 
tions  because  these  papers  assume  that  channel  feedback  information  is  very 
limited,  which  is  characteristic  of  spread-spectrum  systems.  Indeed,  in 
the  appendices  of  this  report  the  traffic  intensity  vs.  packet  error  prob¬ 
ability  tradeoff  (Appendix  B)  and  a  possible  implementation  of  recursive 
retransmission  procedures  as  in  [5]  (Appendix  C)  are  each  given  in  the 
context  of  a  FH-system  such  as  the  Navy's  intra-task-force  communications 
network. 

In  Appendix  D  some  results  from  [9]  are  summarized.  In  this  paper 
the  delay/throughput  tradeoff  of  a  random-access  system  is  studied  under 
the  assumption  of  a  very  limited  amount  of  feedback.  It  is  found  that 
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there  is  a  potential  for  instabilities  if  the  feedback  information  is 
insufficient. 

In  [8]  we  developed  a  numerical  method  for  finding  the  invariant 
distribution  for  a  class  of  Markov  processes.  The  method  is  useful  for 
performance  evaluation  of  certain  random  access  strategies,  as  shown  in 
{10].  In  Appendix  E  the  method  of  [8]  is  outlined  and  some  of  the  results 
from  [10]  are  summarized. 
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I .  INTRODUCTION 

Several  communications  systems  currently  being  developed  have  the 
following  common  features.  Frequency-hopped  spread-spectrum  modulation  is 
employed  with  a  hopping  rate  not  greater  than  the  data  rate.  Multiple-access 
capability  is  required,  because  with  high  probability  two  or  more  terminals 
will  be  transmitting  simultaneously.  During  transmission  the  spread-spectrum 
signals  encounter  severe  fading,  which  causes  reduced  signal  strength  and  may 
produce  intersymbol  interference  or  other  dispersive  effects.  These  systems 
are  described  in  current  terminology  as  slow-frequency-hopped  (SFH)  spread- 
spectrum  multiple-access  (SSMA)  communications  systems  with  fading  channels. 

In  this  paper  we  present  bounds  and  approximations  for  the  average 
probability  of  error  for  SFH/SSMA  communications  over  fading  channels.  Two 
Important  classes  of  fading  models  are  considered:  the  class  of  nonselective 
Rician  fading  channels — which  Includes  the  additive  white  Gaussian  noise 
channel  and  the  nonselective  Rayleigh  fading  channel  as  special  cases--and  the 
selective  wide-sense-stationary  uncorrelated-scattering  fading  channel.  The 
data  modulation  is  binary  frequency-shift  keying  (FSK) ,  but  many  of  the  results 
apply  to  differential  phase-shift-keying  (DPSK)  as  well.  Noncoherent  demodula¬ 
tion  of  the  data  is  employed,  partly  because  we  do  not  require  coherent  frequency 
hopping  and  dehopping.  The  communications  network  is  assumed  to  be  asynchronous; 
that  is,  a  given  terminal  makes  no  attempt  to  coordinate  its  transmissions  with 
those  of  other  terminals.  This  may  be  due  to  the  lack  of  an  accurate  timing 
reference  or  because  of  the  variation  in  propagation  times  among  the  different 
communication  paths  in  the  network.  The  point  here  is  that  even  if  the  trans¬ 
mitters  have  a  common  clock  they  cannot  adjust  their  transmission  times  to 
provide  coordinated  arrival  times  at  all  of  the  receivers  in  the  network. 
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In  the  analysis  of  SFH/SSMA  systems  there  are  two  approaches  to  the 
modeling  of  the  frequency  hopping  patterns:  general  random-process  models  may 
be  employed  or  specific  (deterministic)  sets  of  hopping  patterns  may  be  con¬ 
sidered.  The  random-process  models  are  often  used  in  an  attempt  to  match 
certain  characteristics  of  extremely  complex  hopping  patterns  which  have  very 
long  periods.  Also  random-process  models  serve  as  substitutes  for  deterministic 
models  when  the  comminications  engineer  is  given  little  or  no  information  about 
the  structure  of  the  hopping  patterns  to  be  used  in  the  system.  Both  random 
patterns  and  a  special  class  of  deterministic  patterns  (based  on  Reed-Solomon 
codes)  are  considered  in  this  paper. 

The  results  obtained  in  this  paper  are  bounds  and  approximations  for  the 
bit  error  probability.  These  results  are  useful  for  both  uncoded  FH/SSMA 
systems  and  fully-interleaved  coded  FH/SSMA  systems.  For  coded  systems  which 
employ  random-error-correcting  codes,  full  interleaving  is  usually  necessary  for 
satisfactory  performance.  We  have  also  obtained  results  (similar  to  those  pre¬ 
sented  in  Section  III)  on  the  probability  of  error  for  FH/SSMA  systems  which 
employ  certain  burst-error-correcting  codes  and  "partial  interleaving",  but  this 
topic  is  beyond  the  scope  of  the  present  paper. 

A  brief  outline  of  the  paper  is  as  follows.  The  model  for  the  SFH/SSMA 
system  is  presented  in  Section  II  where  our  models  for  the  various  subsystems 
and  signals  are  described.  The  effect  of  nonselective  fading  on  the  probability 
of  error  in  a  SFH/SSMA  system  is  considered  in  Section  III.  A  more  precise 
analysis  is  given  in  Section  IV  for  the  special  case  in  which  the  channel 
exhibits  nonselective  Rayleigh  fading.  Finally,  selective  fading  is  considered 


in  Section  V. 
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II.  SYSTEM  MODEL 

The  transmitter  for  the  slow- frequency-hopped  spread-spectrum  signal 
is  shown  in  Figure  1.  There  are  K  such  transmitters  in  the  spread-spectrum 
multiple-access  system.  The  k-th  data  signal  bk(t)  is  a  sequence  of 
positive  and  negative  rectangular  pulses  of  duration  T.  The  amplitude  of 
the  4-th  pulse  for  the  k-th  signal  is  denoted  by  (i.e.,  b.  (t)»b^k^ 

for  4T  as  t  <  (4+l)T),  and  b^  is  either  +1  or  -1  for  each  k  and  4.  The 
data  signal  (t)  is  the  input  to  an  FSK  modulator,  and  the  corresponding 
output  is 

ck(t)  -  cos{2TT[fc  +  bk(t)A]t  +  0k(t)}  (1) 

where  A  is  one-half  the  spacing  between  the  two  FSK  tones.  The  signal 

0k(t)  is  the  phase  signal  introduced  by  the  FSK  modulator;  that  is,  if 

bfk^  *  m  then  0,  (t)  ®  0.  for  4T  s£  t  <  (4+l)T  where  0,  is  the  phase  of 
x  K  K,m  R,m 

the  tone  at  frequency  f£  +  mA  for  ■  ■  +1  or  m  «  -1. 

The  FSK  signal  is  then  frequency-hopped  according  to  the  k-th  hopping 
pattern  fk(t)  which  is  derived  from  a  sequence  (f^)  *  . . .  ,f^  ,fQ^,f  ^ ,  ... 
according  to 

fk(t>  -  fj(k),  jlh  S  t  <  (j+l)Th.  (2) 

The  parameter  T^  is  the  time  between  hops  (also  called  the  dwell  time). 

For  slow- frequency -hopping  T^  is  an  integer  multiple  of  T.  The  frequencies 
(k) 

f,  ’  are  all  from  the  set  S  *  fv  :  1  s  n  £  q]  which  is  ordered  such  that 
j  n  n 

vn  <  for  each  n.  Let  A'  be  the  minimum  spacing  between  the  frequencies 
in  the  set  S,  end  let  »  T^/T  be  the  number  of  data  bits  per  hop. 

The  band-pass  filter  shown  in  Figure  1  removes  unwanted  frequency 
components  present  at  the  output  of  the  multiplier.  The  signal  at  the 
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output  of  the  filter  is 


where 


sk(t)  -  V2P  cos[2irfk(t)t  +  cpk(t)]. 


(3) 


fi,(t)  *  f  +  b.  (t)A  +  f.  (t) 


(*) 


and 


9k(c)  “  9k^C)  +  “k^-*' 


(5) 


The  signal  c*k(t)  represents  the  phase  shifts  introduced  by  the  frequency 
hopper  as  it  switches  from  one  frequency  to  another.  Accordingly,  c,  (t) 

/Jr  \ 

is  constant  on  the  time  intervals  that  f.  (t)  is  constant.  Let  or* 
denote  the  value  of  ark(t)  on  the  interval  [jT^, (j+l)T^). 

The  quantity  P  that  appears  in  (3)  is  the  power  of  the  k-th  signal 
at  the  receiver  in  the  absence  of  fading.  In  order  to  account  for  fading, 
we  will  multiply  P  by  a  suitable  factor  to  obtain  the  average  power  in 
the  received  signal.  For  simplicity  we  have  assumed  that  the  signals 
sk(t)  all  have  the  same  power.  However,  as  we  will  point  out  later,  we 
obtain  error  probability  bounds  that  are  valid  even  if  the  power  levels 
are  not  equal. 

Since  we  are  considering  an  asynchronous  multiple-access  system,  we 
allow  an  arbitrary  time  delay  Tk  for  the  k-th  communication  link  (l£k£K). 
Thus  the  received  signals  are  sk(t-Tk),  1  £  k  <  K.  For  the  random  hopping 
patterns  that  will  be  considered  in  subsequent  sections,  it  is  sufficient 
to  consider  time  delays  modulo  T^.  In  order  to  allow  for  the  possibility 
of  deterministic  periodic  hopping  patterns,  we  consider  time  delays  modulo 
NT^  where  N  is  the  period  of  the  patterns  for  deterministic  hopping 
patterns  or  N *  1  for  random  hopping  patterns.  Thus  we  may  restrict 
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attention  to  time  delays  in  the  range  0  £  <  NT^.  Similarly  we  are 

only  concerned  with  phase  angles  modulo  2rr,  so  we  may  restrict  attention 
to  phase  angles  in  the  interval  [0,2rr]« 

The  analysis  presented  in  this  paper  does  not  account  for  adjacent 
channel  interference  in  the  frequency-hopping  system  or  for  interference 
between  the  two  FSK  tones  of  a  given  signal.  Instead  we  are  primarily 
concerned  with  multiple-access  interference  and  the  effects  of  fading  such 
as  intersymbol  interference  and  reduced  signal  strength.  In  order  to 
focus  on  multiple-access  interference  and  fading,  we  made  certain  simpli¬ 
fying  assumptions  concerning  the  frequency  spacings  A  and  A'.  It  is  enough 
for  our  purposes  to  have 

A'  »  A  +  T-1  (6a) 

and 

A  »  T*1.  (6b) 

However,  it  is  possible  to  relax  these  conditions  somewhat,  expecially  for 
nonselective  fading.  For  example  if  the  fading  is  nonselective  then  it  is 
sufficient  to  replace  the  constraint  A  »  T  *  by  the  condition  A  -  m/2T 
positive  integer  m  (the  case  m  ■  1  is  of  greatest  interest).  In  the 
absence  of  time-selective  fading  our  results  are  valid  if  A'  is  about 
3 (A  +  T  S  or  larger,  and  they  are  likely  to  be  fairly  good  approximations 
even  if  A*  «  2 (A  +  T  *).  However,  frequency  dispersion  can  expand  the 
signal  bandwldths  so  that  A'  »  A  +  T  is  needed  for  time-selective 
fading . 

Under  our  assumptions,  the  frequency  band  that  contains  the  signals 
s^t)  is  approximately  the  band  from  f£  +  -  A  to  f£  +  +  A.  The 

center  of  this  band  is  at  frequency  f^  -  fc  +  $(v^-v^).  The  (one-sided) 


bandwidth  W  is  approximately  +  2A.  Under  our  assumptions 

W  »  v  -  v,  i  (q-l)A' . 

q  1 

In  the  absence  of  fading  and  noise  the  received  signal  is  given  by 


K 

s(t)  -  Z  sk(t  -  Tk).  (7) 

k=l 

We  focus  our  attention  on  the  receiver  for  the  i-th  signal,  and  in  doing 
so  we  may  select  the  time  reference  such  that  ■  0.  The  variables  rk 
are  then  delays  (modulo  NT^)  relative  to  this  time  reference. 

The  receiver  for  the  i-th  signal  is  shown  in  Figure  2.  The  received 
signal  s(t),  which  is  a  faded  version  of  s(t),  is  the  input  to  the  first 
band-pass  filter.  This  filter  has  center  frequency  approximately  f^  and 
bandwidth  approximately  W  so  s(t)  is  passed  without  distortion.  This 
filter  is  followed  by  the  i-th  dehopper  which  is  synchronized  in  frequency 
and  time  to  the  i-th  frequency -hopping  signal  f ^ (t) .  The  dehopper 
introduces  a  phase  signal  (3^(t)  which  is  analogous  to  the  phase  signal 
a^(t)  introduced  by  the  frequency  hopper.  The  phase  signal  £L(t)  is 
constant  during  the  time  intervals  between  hops  (i.e.,  when  f^Ct)  is 
constant).  The  constant  value  of  p^(t)  for  jT^  £  t  <  (j+l)Tk  is  denoted  by 


The  time  delays,  phase  angles,  and  data  symbols  are  modeled  as 

mutually  independent  random  variables  each  of  which  is  uniformly  distributed 

on  the  appropriate  set  (cf.  [4]  or  [6]).  The  random  time  delays  are  the 

random  variables  rk*  The  random  phase  angles  that  are  of  primary  interest 

are  9.  ,  or^,  and  An  important  feature  of  our  model  for  asynchronous 

x,m  j  J 

spread-spectrum  multiple-access  systems  is  that  addition  of  phase  angles 
is  modulo-2ir  addition.  This  feature  is  critical  to  our  assertions  concerning 


Figure  2.  Receiver  Model. 
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the  distributions  and  the  statistical  independence  of  the  phase  angles  (the 
basis  for  these  assertions  is  given  on  pp.  159-160  of  [6]). 

The  output  of  the  dehopper  is  then  passed  through  a  band-pass  filter 
which  is  designed  to  remove  certain  unwanted  signals  such  as  the  double¬ 
frequency  components  of  the  i-th  signal  itself,  the  sum  and  difference 
frequency  components  due  to  the  other  K-l  signals  (except,  of  course,  those 
that  happen  to  be  at  the  same  frequency  as  the  i-th  signal),  and  the 
thermal  noise  that  is  outside  the  frequency  band  occupied  by  the  i-th 
signal.  The  bandwidth  B  of  this  band-pass  filter  is  less  than  A'  but 
usually  larger  than  2 (A  +  T  *).  If  (A  +  T  S  «  B  <  A'  then  the  thermal 
noise  present  at  the  output  of  the  band-pass  filter  which  follows  the 
dehopper  has  a  bandwidth  larger  than  that  of  the  FSK  demodulator.  This 
simplifies  the  analysis  of  the  demodulator. 

As  shown  in  Figure  2  the  FSK  demodulator  has  two  branches.  Each  branch 
~2 

forms  a  statistic  where  m=l  corresponds  to  the  upper  branch  and  m=-l 
corresponds  to  the  lower  branch.  Each  of  these  two  branches  has  two 
components.  In  the  in-phase  component  the  signal  is  multiplied  by 
cos [2tt (f^  +  mA)t],  and  the  quadrature  component  it  is  multiplied  by 
sin[2Tr(fc  +  mA)t]. 

Consider  the  reception  of  the  data  bit  The  outputs  of  the 

in-phase  components  of  the  two  branches  are  given  by 

«+l)T 

Zc  m  “  I  r.(t)cos[2TT(f  +rnA)t]dt  (8) 

’  IX 

for  m  *  +  1,  where  r^(t)  is  the  output  of  the  band-pass  filter  which 
follows  the  i-th  dehopper  (i.e.,  r^(t)  is  the  input  to  the  i-th  FSK 
demodulator).  Notice  that  in  general  Z£  m  depends  on  both  l  and  i. 
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However,  if  the  random  hopping  patterns  are  stationary  and  identically 

distributed  and  the  fading  process  is  stationary  and  not  frequency 

selective,  then  the  distribution  of  the  random  variable  Z  will  not  depend 

c,m 

on  either  l  or  i.  In  case  the  hopping  patterns  are  deterministic  or  the 
fading  is  frequency  selective  then  we  provide  upper  bounds  on  the 
probability  of  error  which  are  independent  of  l  and  i.  The  outputs  of 
the  quadrature  components  of  the  two  branches  are  denoted  by  Z  for 

S  ,111 

m  *  +  1.  The  random  variables  Z  ,  which  are  defined  by  (8)  with 
—  s,m*  J 

cos[>]  replaced  by  sin[*],  have  the  same  properties  as  Z 

c  ,m 
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III.  PERFORMANCE  OF  FH/SSMA  SYSTEM  WITH  NONSELECT IVE  FADING 

The  channels  considered  in  this  section  are  the  nonselective  slow- 
fading  channels.  For  the  frequency-hopped  spread-spectrum  system  described 
in  the  previous  section  this  means  that  the  signal  at  the  input  to  the 
first  band-pass  filter  in  the  i-th  receiver  (see  Figure  2)  is 

K 

r (t)  -  n(t)  +  Z  y.  (t-Tk),  (9) 

k-1 

where  for  (T  s  t  <  (jC+1)T  the  signal  yk(t)  is  given  by 

yk (t )  »J2 P  Ajk)  cos  [2trfk (t  )t  +  ^  (t )  +  tjk)  ] .  ( 10 ) 

The  thermal  noise  n(t)  is  white  Gaussian  noise  with  spectral  density 

^-Nq.  Notice  from  comparisons  of  (9)  and  (10)  with  (7)  and  (3),  respectively 

that  yk(t)  is  a  faded  version  of  sk(t)  and , in  the  absence  of  noise »r(t)  is 

s(t)  *  yk ( t  " T k)  which  is  a  faded  version  of  s(t). 

The  amplitude  of  the  fading  signal  yk(t)  during  the  time  interval 

(k) 

iX  i  t  <  Oi+l)T  is  represented  by  a  nonnegative  random  variable  A^  ,  and 

(k) 

the  phase  shift  due  to  the  fading  is  denoted  by  .  In  this  section  the 
only  assumption  that  we  make  concerning  the  signal  amplitudes  is  that  they 
are  constant  during  the  data  bit  interval.  The  sequence  of  amplitudes 
(A^  )  *  . . . ,A_^  >  Aq  ,  A^  ,  ...  may  be  any  stationary  random  sequence. 

In  particular  we  place  no  restrictions  on  the  statistical  dependence  of 
amplitudes  in  different  data  bit  intervals.  Consider  the  set 
{ A^k^  :  jN^  £  l  <  (j+lJN^}  of  amplitudes  for  the  data  bits  that  are 
transmitted  during  the  j-th  hopping  interval  [jT^, (j+l)Th) .  This  interval 
contains  the  data  bits  for  jN^  £  l  <  (j+lJN^.  Among  the  cases  of 
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(k)  (k) 

interest  are  the  two  extreme  cases  (i)  A.  “A  for  all  J l  and  m  in  the 

L  m 

Ck')  fk) 

same  hopping  interval  and  (ii)  A^  and  A^  are  independent  if  l  t  m 

but  i  and  m  are  in  the  same  hopping  interval.  Case  (i)  corresponds  a 

system  with  no  interleaving  and  a  channel  with  slow  fading  relative  to  the 

hopping  rate.  An  example  of  case  (ii)  arises  in  a  system  which  is  fully 

interleaved  (e.g.  if  a  random-error-correcting  code  is  to  be  employed). 

Although  these  are  the  two  specific  cases  of  greatest  interest,  there  is 

no  need  to  restrict  attention  to  such  special  cases  in  this  section. 

(k) 

Similarly,  we  impose  no  restrictions  on  the  phase  sequence  ($£  ') ;  all 
that  is  required  is  a  constant  value  for  the  phase  during  the  data  bit 
intervals.  Notice  from  (l)-(5)  that  for  S  t  <  (X+1)T  the  phase  of  the 
signal  in  (10)  is  given  by 


9,  +  Of. 

k,m  j 


CO 


(11) 


where  j  is  the  integer  part  of  A/N^-  Under  quite  general  conditions  the 
(k) 

phase  0'  '  is  uniformly  distributed  on  [0,2tt]  because  the  addition  in  (11) 

If 

is  modulo-2TT.  For  example,  it  is  enough  to  assume  that  one  of  the  phase 
angles  which  appears  on  the  right-hand  side  of  (11)  is  uniformly  distributed 
and  that  they  are  mutually  independent  (see  pp.  159-160  of  [6]  for  the 
rationale  for  this  statement). 

There  are  two  different  phenomena  which  contributed  to  errors  in  the 
system  under  consideration.  First,  even  in  the  absence  of  noise  and  fading, 
errors  may  occur  in  a  frequency-hopped  spread-spectrum  multiple -access 
system  when  a  signal  is  hopped  to  a  frequency  slot  that  is  occupied  by 
another  signal.  Whenever  two  different  signals  simultaneously  occupy  one 
frequency  slot  we  say  a  hit  occurs.  Second,  even  in  the  absence  of  hits, 
errors  may  occur  due  to  the  fading  and  additive  noise.  The  first  step  in 
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analyzing  the  overall  probability  of  error  ia  to  evaluate  the  probability 
of  a  hit  for  various  types  of  hopping  patterns. 


Probability  of  a  Hit 


Consider  as  before  the  i-th  receiver  during  reception  of  the  X-th 
data  bit.  For  a  nonselect ive  fading  channel  we  say  that  a  hit  from  the 
k-th  signal  occurs  during  the  X-th  data  bit  if 


fk(t“Tk}  "  fi(t)  (12) 

for  at  least  one  value  of  t  in  the  X-th  data  bit  interval  [AT, (X+1)T) . 

As  pointed  out  in  Section  II,  we  can  let  N  »  1  in  considering  stationary 

random  hopping  patterns.  It  follows  that  the  probability  9^  of  a  hit 

from  the  k-th  signal  during  the  X-th  data  bit  interval  does  not  depend 

on  l  for  such  patterns.  If  the  K  hopping  patterns  {  (f^)  :lSk«  K} 

Ck') 

are  also  mutually  independent  and  identically  distributed  then  9:  does 

z 

not  depend  on  k  either,  and  hence  we  denote  it  by  9  for  such  patterns. 

We  first  consider  two  different  models  for  stationary  random  hopping 

patterns  and  give  the  value  of  9  for  each  case. 

(k) 

Suppose  the  random  process  (fj  ')  is  a  stationary  Markov  process  with 
transition  probabilities  given  by 


P(fj+l  ’vnlfj(k)  ’  V  "  <q_irl 


for  l  £  n  £  q,  1  £  r  at  q,  and  n  j*  r.  It  follows  that  for  these  patterns 


PCfjil  -  fj(k))  -  0 


and  hence 
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*  ■  ^  +  • 
:(k). 


(15) 


Because  of  (13)  the  process  (fj  ;)  is  a  random  process  with  first-order 
distribution  given  by 


P(fj(k)  »  VQ)  -  q‘\  1  *  n  *  q. 


(16) 


If  instead  of  (13)  we  consider  random  hopping  patterns  for  which  f 
is  independent  of  f^k^  and  the  distribution  of  fjk^  £s  given  by  (16)  for 
each  j,  then  (14)  should  be  replaced  by 


(k) 

J+l 


P(f^  -  fjk))  =■  q-2 


(17) 


and  thus  the  probability  of  a  hit  is 


(18) 


Notice  that 


9  &  (1  +  ~) 

q  ^ 


(19) 


and  if  q  is  large  then 


q  Njj 


(20) 


(cf.  equation  (15)).  Thus  for  large  q  these  memoryless  hopping  patterns 
give  approximately  the  same  probability  of  a  hit  as  the  first-order  Markov 
patterns . 
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In  general  for  a  set  of  deterministic  hopping  patterns  the  probability 

(k) 

depends  on  both  k  and  j t.  One  set  of  deterministic  hopping  patterns  that 
has  very  good  properties  is  derived  from  a  Reed-Solomon  code,  so  we  refer  to 
it  as  a  set  of  Reed-Solomon  (RS)  patterns  [7].  Given  a  prime  number  q  of 
frequency  slots,  the  particular  set  of  RS  patterns  of  interest  here  consists 
of  N  *  q-1  sequences  of  period  N  (of  course  we  can  always  choose  a  subset  if 
fewer  patterns  are  needed) .  Each  sequence  is  nonrepeating; 
that  is,  for  each  sequence  (f^),  6(fj,fQ)  *  0  for  n  #  j  and  0^  nS  N-l, 
where 


6 (u,v) 


u  *  v 
U  +  V  . 


(21) 


The  property  of  RS  patterns  that  is  of  primary  importance  here  is  that  for 
any  two  patterns  (fj  )  and  (fj1^). 


V  § (f^k) ,f 51))  ^  1 

n*0  n  1 


(22) 


for  each  j.  Property  (22)  is  actually  valid  for  any  set  of  nonrepeating 
patterns . 

Since  is  uniform  on  [0,NTb]  for  k  +  i,  then  it  follows  from  (22)  that 


*00  <*  k  hH  +  a+i-) 

NTh  N  U  +  Nb}  q-1  (1  +  Nb)  * 


(23) 


(k) 

Actually  (22)  implies  the  stronger  statement  that  either  ^  *  0  or 

(k) 

9^  •  9  .  Since  the  number  of  frequency  slots  q  is  larger  than  the 

(k) 

period  N  *  q-1,  then  it  is  possible  to  choose  the  hopping  pattern  (f^  ) 

such  that  *  0  for  Nb  different  values  of  j t  in  the  range  0  <  l  <  NNb 
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Notice  that  for  large  q 

*  a+ir>  w 

q  Nb 

is  a  good  approximation  for  the  upper  bound  (cf.  (15)  and  (20)). 

Of  primary  interest  for  our  subsequent  analysis  is  the  probability 

A 

9^  of  one  or  more  hits  from  the  K-l  signals  (corresponding  to  k  +  i)  during 

A 

the  Z-th  data  bit  interval.  For  stationary  random  patterns  9^  does  not 

A 

depend  on  i  so  we  denote  it  by  &.  If  the  patterns  are  also  mutually 
independent  and  identically  distributed  then 

9  -  1  -  (l-^)*”1  ,  (25) 

where  9  is  the  probability  of  a  hit  from  a  given  signal.  For  the  first-order 
Markov  patterns  (25)  and  (15)  imply 

9  -  1  -  Cl  -  J  (1 +£-)}*’ 1  .  (26) 

b 

If  the  patterns  are  sequences  of  independent  random  variables  (i.e. 
memoryless  patterns)  satisfying  (16)  then 

9  •  1  -  Cl  -  i[l+^-(l  -i)]}K_1  .  (27) 

Next  we  consider  the  probability  ^  of  one  or  more  hits  in  the  i-th 
data  bit  interval  for  deterministic  patterns.  Since  the  random  variables 
k  +  i,  are  mutually  independent,  then  for  any  deterministic  hopping 
patterns 

»  .1.  (  n  U  -*<*>1 
i  }  k-l  1 
lk*i 


(28) 
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For  RS  patterns  (23)  implies 

^  1  -  [1  -  *  ]K_1  “  1  -  U  -  (l  +  ^-)}K'1  .  (29) 

b 

A 

where  the  symbol  9  is  used  to  denote  an  upper  bound.  Notice  from  (27)  and 
(29)  that  for  large  q 

1  -  {1  -  Jd  +  J-)}**1  (30) 

q  «b 


for  the  sequences  of  independent  random  elements  and  the  RS  sequences. 

Notice  from  (26)  that  the  expression  given  in  (30)  is  the  exact  value  of  9 
for  the  first-order  Markov  patterns. 

B .  Bounds  and  Approximations  for  the  Probability  of  Error 

For  a  nonselective  fading  channel  the  bit  error  probability  P  in 

&  yZ 

a  slow- frequency-hopped  spread-spectrum  multiple-access  comuunications  system 
can  be  written  as 


e,2 


V1  -*i>  +  h,i 


(31) 


where  Pq  is  the  conditional  probability  of  error  for  the  2-th  bit  given  that 
there  are  no  hits  and  P^  ^  is  the  conditional  probability  of  error  for  the 
2-th  data  bit  given  there  is  at  least  one  hit.  Notice  that  Pq  does  not 
depend  on  2.  In  general  P^  ^  depends  on  2  but,  as  will  be  seen  from  the 
numerical  results,  it  is  sufficient  for  many  purposes  to  use  the  bounds 


0  <  P 


1.1 


<  fc. 


A 

Recall  that  for  stationary  random  hopping  patterns  9^  does  not  depend 

A  A 

on  i  (and  hence  it  is  denoted  by  9) .  For  RS  patterns  9^  depends  on  i  but  its 

A 

upper  bound  9  given  by  (29)  does  not.  Hence  for  all  of  these  patterns  we 
have  the  lower  bound 

Pe,lSP L'V1-^  <32> 

and  the  uppeT  bound 

Pe,i£  Pu  S  h  +  ‘  P0+  <*-po*  ‘33> 

where  P  is  given  by  (26),  (27),  or  (29),  depending  on  which  type  of  hopping 
patterns  are  employed.  The  lower  bound  is  the  same  as  we  previously  presented 
in  [5],  but  the  upper  bound  of  (33)  is  a  slight  improvement  of  the  upper 
bound  presented  in  [5]. 

It  is  tempting  to  use  P  .  2  P.  in  place  of  (32),  and  we  certainly 

fij*  '-I 

believe  this  tighter  lower  bound  to  be  valid  for  independent  time  delays, 
data  streams,  and  hopping  patterns.  Under  these  conditions  it  is 
intuitively  clear  that  multiple-access  interference  cannot  decrease  the 
average  probability  of  error.  However,  the  lower  bound  of  (32)  has  the 
advantage  that  it  holds  under  more  general  conditions  (such  as  for  dependent 
time  delays,  data  streams ,  and  hopping  patterns). 

The  bounds  given  in  (32)  and  (33)  are  valid  even  if  the  power  levels  are 
not  the  same  for  the  various  signals  or  the  hopping  patterns  are  statistically 
dependent.  As  might  be  expected,  the  imposition  of  additional  restrictions 
on  the  system  leads  to  more  precise  results.  In  Section  IV  we  present  such 
results  for  a  more  restrictive  channel  model.  However,  even  with  the  full 
generality  of  the  nonselective  fading  channel  model  considered  in  this 
section,  we  can  improve  the  lower  bound  and  obtain  a  useful  approximation 
if  we  consider  equal  power  signals  and  add  certain  constraints  on  the 
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hopping  patterns  and  the  binary  data  streams.  The  hopping  patterns  are 
assumed  to  be  stationary,  nutually  Independent,  Identically  distributed  random 

patterns,  and  the  data  sequences  are  stationary,  memory less.  Independent 

GO 

random  sequences  with  distribution  given  by  P(b^  =  m)  =  %  for  m  =  +l  and  m*-l. 

The  lower  bound  can  be  improved  for  such  systems  by  providing  a  nonzero 

A 

lower  bound  for  the  term  P^  ^9^  of  (31).  One  such  bound  is  obtained  as  follows. 
Consider  the  conditional  probability  of  error  in  the  X-th  data  bit  given  a 
"full"  hit  from  the  k-th  signal  (i.e.,  given  that  (12)  holds  for  all  t  in 
(XT,(X+1)T))  and  given  the  k-th  signal  transmits  -b^  for  the  two  consecutive 
bit  intervals  of  interest.  This  conditional  probability  of  error  is  equal  to 
The  conditional  probability  of  a  "full"  hit  (given  a  hit  has  occurred)  is  not 
smaller  than  (Nb-1) / (Nfe +  1) ,  and  the  probability  of  two  consecutive  transmissions 
of  a  particular  tone  is  Finally,  we  use  the  fact  that  (25)  implies 

9l  *  (K-l)*»(l  -<?)K'2, 


which  is  just  the  statement  that  the  probability  of  one  or  more  hits  is  not  less 
than  the  probability  of  exactly  one  hit.  From  the  above  we  conclude  that 


(Nb-D 

p  &  i  — - - 

l,x  X  a(Nb+i) 


(K-l)  9  (1  - 9 ) 


K-2 


so  that  the  improved  lower  bound  is 


PL  + 


(N  -1) 

4-ri)(K-i)*(i  -*> 


K-2 


(34) 


tfe  use  tilde  ( ~  )  to  denote  bounds  and  approximations  which  are  valid  for  the 
restricted  class  of  systems  only  (i.e.,  equal  power  signals,  memoryless 
independent  data  sequences,  independent  hopping  patterns). 


An  approximation  which  is  valid  under  the  same  conditions  is 
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?e,l  ~  \  ^  PL  +  h  ^+P  0)(K_1)^(1  -^)K'2  •  (35) 

This  approximation  is  very  accurate  whenever  q/K  is  large  because  it  is 
based  on  the  assumption  that  the  probability  of  a  multiple  hit  (i.e.  hits 
from  two  or  more  signals  in  a  given  data  bit  interval)  is  negligibly  small 
in  comparison  to  the  probability  of  a  hit  from  only  one  signal. 

Comparisons  of  the  bounds  and  the  approximation  are  given  in  Table  1  for 
various  values  of  Pq,  K,  q,  and  N^.  The  hopping  patterns  are  the  first-order 
Markov  patterns  for  the  data  in  Table  1,  but  in  view  of  (30)  the  results 
would  not  be  significantly  different  for  the  other  patterns  described  above. 

C .  The  Nonselectlve  Rician  Fading  Channel 

The  bounds  and  approximation  given  in  (32) -(35)  can  be 
applied  to  any  particular  nonselectlve  fading  channel  by  substituting  the 
appropriate  expression  for  PQ  in  these  results.  In  this  section  we  consider 
the  Rician  nonselectlve  fading  model  in  which  each  transmitted  signal  results 
in  a  received  signal  that  is  the  sum  of  a  nonfaded  version  of  the  trans¬ 
mitted  signal  and  a  (nonselectlve)  Rayleigh  faded  version  of  the  transmitted 
signal.  The  difference  in  the  propagation  times  for  these  two  components 
is  sufficiently  small  compared  with  the  data  bit  duration  T  that  the  overall 
channel  is  nonselectlve.  This  model  is  discussed  in  [9]  where  the  nonfaded 
component  is  called  the  fixed  or  specular  component  and  the  Rayleigh-faded 
component  is  called  the  random  or  scatter  component.  In  some  applications 
the  nonfaded  component  arises  from  a  direct  path  between  the  transmitter 
and  the  faded  component  arises  fr.om  a  reflection. 


'  I 


Table  1.  Lower  bounds,  approximation,  and  upper  bound 
on  the  probability  of  error  for  a  FH/SSMA 
system. 


a)  K  «  15,  q  *  1000,  and  »  10 


po 

PL 

PL 

PA 

pu 

0.100 

0.098 

0.100 

0.103 

0.106 

0.050 

0.049 

0.051 

0.053 

0.057 

0.030 

0.030 

0.031 

0.034 

0.037 

0.020 

0.020 

0.021 

0.024 

0.027 

0.010 

0.010 

0.011 

0.014 

0.017 

0.005 

0.005 

0.006 

0.009 

0.013 

K  =  15, 

P0 

q  *  100, 

PL 

and  *  5 

PL 

PA 

pu 

0.100 

0.084 

0.096 

0.128 

0.162 

0.050 

0.042 

0.054 

0.081 

0.120 

0.030 

0.025 

0.037 

0.063 

0.103 

0.020 

0.017 

0.029 

0.054 

0.095 

0.010 

0.008 

0.020 

0.045 

0.086 

0.005 

0.004 

0.016 

0.041 

0.082 

K  -  25, 

po 

q  *  250, 

PL 

and  ■ 

PL 

hO 

o 

> 

PA 

0.100 

0.09Q 

0.100 

0.118 

0.138 

0.050 

0.045 

0.056 

0.070 

0.093 

0.030 

0.027 

0.037 

0.051 

0.075 

0.020 

0.018 

0.028 

0.042 

0.066 

0.010 

0.009 

0.019 

0.032 

0.057 

0.005 

0.005 

0.015 

0.027 

0.053 

i 
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The  amplitude  S  of  the  sum  of  the  two  components  of  the  received  signal 

is  a  random  variable  with  a  Rician  distribution  (see  [8]  or  [9]).  Since  we 

are  interested  in  the  conditional  probability  of  error  given  there  are  no 

hits,  we  can  assume  in  all  chat  follows  that  only  the  components  of  the 

i-th  signal  are  present  at  the  i-th  receiver  (during  the  data  bit  interval 

under  consideration).  Let  p  be  the  normalized  bit  energy  to  noise  density 
2 

ratio,  so  that  S  p  is  the  actual  received  energy  to  noise  density  ratio. 

2 

Hence  for  noncoherent  FSK  the  probability  of  error  given  S  *  a  is  \  exp(-%a  p). 

For  the  Rician  channel  the  density  function  f  for  the  amplitude  S  is 

b 

fg(a)  »  (a/a2)exp{-%(a2 +a'2)/o2}l0(Qf  a/CT2)  (36) 

2 

for  a  >  0,  where  a  represents  the  strength  of  the  nonfaded  component, 

2 

2a  is  the  expected  value  of  the  strength  of  the  faded  component,  and  Iq 
is  the  zero-th  order  modified  Bessel  function.  The  average  probability  of 

ic 

error  for  noncoherent  FSK  is  [9] 


^Corresponding  results  for  binary  DPSK  are  obtained  by  replacing  p 
by  2p  in  (37) . 
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Then  we  can  write 


pn  =  ~2— ^-'2^—  exPt"A/  [y2A  +  2(y2  +  1)  ]}  .  (39) 

yA  +  2(y  +1) 


2  2  2  2 

Two  limiting  cases  of  interest  are  a  =0  and  a  =0.  If  a  =  0  (y  =0) 

then  there  is  no  faded  component,  and  the  channel  is  just  an  additive  white 

2 

Gaussian  noise  channel.  In  this  case  A  =  a  p,  and  the  probability  of  error 
reduces  to 

PQ  =  k  expC-%  A}  .  (40) 

2  2 
If  a  =0  the  channel  is  a  nonselective  Rayleigh  fading  channel,  A  =  2a  p , 

and  the  probability  of  error  is 


P  —  "  '  ~ 1 
0  A  +  2 


(41) 


2  2 

An  examination  of  (39)  as  a  function  of  y  shows  that  for  y  =  10  the 

probability  of  error  for  Rician  fading  is  nearly  the  same  as  for  Rayleigh 

-2  2 

fading.  For  example,  if  A  is  12  dB  then  Pq  is  1.81  x  10  for  y  =*  0, 

4.41  x  10"2  for  y2  =■  0.1,  4.53  x  10'2  for  y2  *  1.0,  and  5.58  x  10-2  for 

2  2-2 
y  =  10.0.  The  value  of  Pq  for  Rayleigh  fading  (y  =  •)  is  5.60  x  10 

In  order  to  apply  (36)-(41)  to  the  slow-frequency-hopped  spread -spectrum 

multiple-access  system,  consider  first  the  expressions  (9)  and  (10)  for  the 

(VO 

received  signal.  The  amplitudes  A)  are  random  variables  with  a  density 

J ft 

function  of  the  form  given  in  (36).  In  general  the  parameters  or  and  a  may 
depend  on  i,  in  which  case  A  and  y  also  depend  on  i.  The  probability  Pq 
Chen  depends  on  i  and  is  given  by  (37)  with  a  and  a  replaced  by  and 
or  by  (39)  with  y  and  A  replaced  by  y  and  It  follows  from  (9)  and  (10) 

that  the  parameter  p  is  given  by  p  ■  PT/Nq. 


* 
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The  next  step  Is  to  substitute  for  Pq  in  (32)-(35)  using  the  expressions 
(37)  or  (39).  If  Pq  depends  on  i  the  bounds  of  (32)  and  (33)  are  valid,  but 
of  course  they  will  also  depend  on  i.  Notice  that  if  a  and  a  depend  on  i, 
then  the  average  power  in  the  received  signal  also  depends  on  i.  That  is, 
the  signals  are  not  required  to  have  equal  power.  The  approximation  given 
in  (35)  is  also  valid  even  if  Pq  depends  on  i,  provided  that  and 

a  for  all  k. 

In  Figure  3  the  approximation  P  ,  which  is  given  by  (35)  with  Pq  replaced 

by  the  expression  in  (39),  is  shown  as  a  function  of  A  =  5/Nq  for  various  values 
2 

of  y  .  For  the  data  presented  in  Figure  3,  the  values  of  a  and  a  (and 
hence  y  and  A)  do  not  depend  on  i.  Additional  numerical  data  can  be  obtained 
from  Table  1  by  evaluating  Pq  from  (37)  or  (39).  Notice  that  for  Rayleigh 
fading  with  <?/Nq  less  than  20  dB  the  value  of  Pq  is  less  than  0.01.  From  Table  1 
we  see  that  for  Pq  ^  0.01,  the  value  of  P^  is  always  less  than  2Pa  and  the  value 
of  is  always  less  than  2PL  for  the  values  of  K,  q,  and  N^  considered  in 
Table  1.  For  K  =  15,  q  *  1000,  and  N^  =  10  we  see  that  for  Pq  —  0.01,  we 
always  have  P  ^  1.2  P  and  P  <  1.25  P  .  Thus,  for  Rayleigh  fading  or 

A  la  U  A 

2 

Rician  fading  with  Y  2  1,  the  bounds  and  approximations  given  in  this  section 
are  sufficiently  accurate  for  the  design  of  slow- frequency-hopped  spread- 
spectrum  multiple-access  systems.  Further  evidence  of  this  is  given  in  the 


next  section. 


Bit  error  probability  (P^)  for  FH/SSMA  system  with  Rician  fading 
(K  =  15,  q  =  1000,  and  N  =  10). 
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IV.  NONSELECTIVE  RAYLEIGH  FADING 

In  this  section  we  present  a  more  exact  analysis  of  the  effects  of 

multiple-access  Interference  and  nonselectivre  fading  for  the  special  case 

in  which  the  fading  is  Rayleigh.  This  analysis  provides  a  more  accurate 

approximation  and  a  tighter  upper  bound  for  the  probability  of  error  than 

is  obtained  by  specializing  the  results  of  Section  III  to  Rayleigh  fading. 

The  system  and  channel  models  are  as  presented  in  Section  III,  and  the 

received  signal  is  as  given  in  (9)  and  (10). 

Since  we  are  considering  only  Rayleigh  fading  in  the  present  section,  the 
(k) 

random  amplitudes  A,  ’  have  a  Rayleigh  distribution.  The  density  function 
(k) 

for  A^  is  given  by  (36)  with  a  *  0  and  a  -  a^.  general  the  second 

moments  =  2 o^  are  different  for  different  signals.  For  the  analysis 

presented  in  this  section  we  assume  that  the  fading  for  different  signals 

is  statistically  independent.  Stated  precisely,  the  requirement  is  that 

A^  ,A^  , . . .  ,a£K^  are  mutually  independent  for  any  choice  of  t  ,l  ,  . . .  ,l  . 

*2  ZK  1  1  K 

The  starting  point  for  the  analysis  of  the  receiver  is  (8).  Since  in 

practice  »  T  for  a  spread -spectrum  system,  the  high  frequency  terms 

in  the  integrand  of  (8)  may  be  ignored.  The  output  of  the  Integration  at 

the  sampling  instant  is  then  given  by 


Z 


c,m 


D  +  (P/8)^  T  £  I^’1*  +  N 

c,m  c,m  c,m 


(42) 


The  first  term  D  is  the  component  due  to  the  signal  s, (t).  If  the 
c  ,m  i 

transmitted  data  bit  is  b^^  for  X  ■  jN^  +  p  then 

Dc,m  "  (P/8)%TAXi)6(bXi),m)c0s[ei,m  +  aj(i)  ‘Pj(1) +4Xi)l  *  (43) 

Since  the  component  is  the  output  of  the  integrator  in  the  absence  of 
multiple-access  and  channel  noise,  it  is  called  the  desired  signal  component. 


(k  i> 

The  multiple-access  interference  lv  from  the  k-th  signal  depends 

c ,  m 

upon  the  delay  t^.  For  convenience  let  Ak  -  and  n^  «  |_(Tk  "  AkTh)/Tj  * 

(k  i) 

where  [yj  denotes  the  integer  part  of  the  real  number  u.  Then  1^  ^  can 
be  expressed  as 

Ic!;1)  =  d<V[^(^)VVVC°8*’  ^k’^  +AL(nk)e2(1’k,nk)cOS  ^k’V1 

(44) 

for  0<  n^  <  p.  The  following  expressions  define  the  various  terms  in  (44). 
First  we  have 

d(A)  -  fi (fj-i*fjt>>  (45) 

for  0^  i  <  N.  Second,  if  L(n)  *  (j  -i^)N^  +  p-n  then 

e^jt.n)  *  6(b^1),m)[Tk-ATh-nT]/T.  (46) 

and 

e2U,n)  -  6(b^),m)((tH-l)T-Tk  +  ^Th]/T  .  (47) 

Finally,  if  b™^  -  m*  and  b^}  *  m"  then 

♦’(*•“)-  ®k. -  •  +  *  P  j(1>  -  *  «  I V’  *  +  f ! $ ]rk  +  *<&.  L)  («) 

and  ijr"(A,n)  is  given  by  (48),  with  m1  replaced  by  m"  and  L(rtf-l)  replaced 
by  L(n).  For  p  <  equation  (44)  is  replaced  by 

l?,«)  •  d<1k<'1)|ALcik+l)'l(,k'“k)c°**  '  < V1- V  +  '2  <  VV“»’  "  'V1 ' 

(49) 

The  only  remaining  case  is  ■  p  for  which  we  have 

1 c^,1'  •  '  <Vl'">  +  d<VA£<J),Z«k-',>c0,*"<V>,)  ■ 

(50) 
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(k)  (k) 

Notice  that  if  we  set  A£^  *  1  and  ^L(n)  *  ®  in  the  above  expressions, 
then  we  obtain  expressions  for  the  in-phase  components  of  the  desired  signal 
and  the  multiple-access  Interference  for  a  system  with  an  additive  white 
Gaussian  noise  channel. 

The  remaining  component  of  Z£  m  is  the  component  Nc  m  which  is  due  to 

the  channel  noise  process  n(t).  It  is  easy  to  show  that  N  is  a  zero-mean 

c ,  m 

Gaussian  random  variable  with  variance  NqT/16. 

The  quadrature  components  are  defined  by  expressions  which  are  analogous 

to  (42)-(44).  In  fact  Z  and  N  are  defined  in  the  same  way  as  above,  and 

s,m  s,m  J 

(k  i) 

the  only  change  that  must  be  made  in  the  definitions  of  D  and  I  is 

J  s,m  s,m 

that  cos(*)  should  be  replaced  by  -sin(*)  in  (43),  (44),  (49)  and  (50). 

We  next  consider  the  average  probability  of  error  where  the  average  is 
computed  with  respect  to  the  phase  angles,  time  delays,  and  data  symbols. 


We  start  by  assuming  that  the  transmitted  data  bit  is  b 


(i) 


+1  where 


X  =  jN^+p  as  before.  Also  the  probabilities  and  expectations  below  are  all 

) 

conditioned  upon  the  data  sequences  (b^  /)  and  time  delays  t^.  For  m  *  +1  or  -1, 
2  2 

let  a  and  a  be  the  variances  of  the  in-phase  and  quadrature  components 
c,m  s,m  ^ 

Zc  m  and  Zg  m  respectively.  Since,  as  we  discuss  below,  these  components  are 
Gaussian  random  variables  with  equal  variances,  the  probability  of  error  is 
given  [8,  p.  587]  by 


2  2  2 
<-!<i  +  <-l> 


(51) 


for  slow  nonselective  Rayleigh  fading  and  noncoherent  FSK  detection.  Under 

the  assumptions  about  the  fading  model  that  were  made  above,  the  desired 

signal  component  D  .  is  a  zero-mean  Gaussian  random  variable  with  variance 
2 

(PT  /16)n. .  Also  notice  that  D  .  m  0. 

X  C  j  •  1 


o. 
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In  order  to  proceed  further  In  the  analysis  of  the  multiple-access 

Interference,  we  need  to  consider  the  nature  of  the  statistical  dependence 

between  aP6^  and  aP*)  and  between  and  for  2  and  2+ 1  In  the  same 

hopping  interval.  These  are  the  random  variables  which  describe  the  fading 

during  adjacent  data  bits.  We  consider  the  two  extreme  cases  described  in 

HO  (k) 

Section  III:  (i)  the  fading  is  constant  in  the  sense  that  A^  *  ^2+1  and 
(lO  (k) 

whenever  2  and  2+1  are  in  the  same  hopping  interval  and  (ii)  the 

fading  is  independent  for  adjacent  data  bits  in  the  same  hopping  interval. 

Under  our  assumptions,  the  multiple-access  interference  component 

l(k,i)  a  Zero-mean  Gaussian  random  variable  with  variance  %  p.,a2(k,i).  For 
c ,  in  Km 

constant  fading,  as  described  by  case  (1)  above,  we  have 

<4<k,i)  -  d^Me^J^n^+e^.r^)]2  (52) 

for  0  <  n^  <  p, 

a2(k,i)  -  d(2k+l)[e1(2k,nk)]2  +  d(2k) te2(2k,nk) ]2  (53) 

for  =  p,  and 

a2(k,i)  -  d(2k+l)(e1(2k,tik)  +  e2(2k>nk)]2  (54) 

for  p  <  nk  <  .  For  independent  fading,  as  described  by  case  (ii)  above, 

we  have 

<j2(k,i)  -  d(2k)C[e1(2k,nk)]2+[e2(2k,nk)]23  (55) 

2 

for  OS  <  p.  If  nk  ■  p,  <7n(k,i)  is  given  by  (53),  and  for  p  <  nk  < 

o2(k,i)  -  d(2k+l)([e1(2k,nk)l2  +  [•2Uk,t»k)  ]23  •  (56) 
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Because  of  the  independence  of  the  fading  for  different  signals,  the  random 
fk  i) 

variables  I  are  independent  when  conditioned  on  the  data  bits  and  time 

c,m 

delays.  As  a  result 


m  »  (PT2/16)t6(l,m)n  +  2  ^.a2(k,i)]  +  NnT/16  . 

c,ra  1  1cm  U 

By  symmetry  we  see  that  a  *  a  .  Thus  (51)  can  be  written  as 
J  J  J  s,m  c,m 


(57) 


<,/V'1  +  k51  W 


l  +  2(  tf/lO"1*  £ 

0  kH 


M'kM’i1[CT-i(k,i)+ai(k»i)1 


(58) 


where  <3  -  f^PT  is  the  energy  per  bit  for  the  received  signal  (in  the  absence 
of  nultiple-access  interference). 

In  order  to  evaluate  the  average  probability  of  error  P^,  we  must  average 
the  expression  in  (58)  with  respect  to  the  time  delays  and  data  symbols.  This 
is  of  course  a  difficult  computation  since  it  involves  the  evaluation  of  K-l 
dimensional  integrals.  However  we  can  obtain  an  approximation  P^  and  an  upper 
bound  Py  which  are  relatively  easy  to  compute.  This  Is  accomplished  by 
observing  that  P£  depends  on  only  through  t^,  and  n^  where 
t^  *  (i^  -  jfc^Th  -  tijcT)/T.  We  can  thus  obtain  a  discrete  approximation  to  the 
integral  with  respect  to  t^  by  approximating  the  uniform  distribution  on  [0,1] 
by  the  discrete  distribution  with  probability  mass  j”2  at  points 

J*^,2J  \ _ , (j-l)j’2  and  probability  mass  (2J)’2  at  the  end  points  0  and  1. 

We  find  that  for  the  first-order  Markov  patterns  and  constant  fading  (case  (i)) 


Pfam(k,i)  -  jj”1}  -  Pj,  0<j<J, 


(59) 
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where  Che  quantities  p  are  defined  completely  by  the  fact  that  their  sum 
is  1  and 


(2Jq)"1(l  +  Nb“1), 


j  =  1,2, .. . , J-l, 


(4Jq)”^(l  +  Nf1)  +  (4q)'1(l-N;1),  j  *  J  . 


For  independent  fading  (case  (ii)) 


PC°m(k.i) 

*  “  Pj  .  °-  J  - 

p{>m(k,i) 

-  [j2+(j-j)2iV1}  -  qj  , 

where  ; 

Pj  and  q^  are  defined  by 

(  (2Jq)“1(l  +  N"1), 

j  -  1,2, .. . ,J-1  , 

v  ■ 

i 

‘  (2Jq)-1  , 

J  -J, 

1  (2jq)”1(l-Nb1), 

j  »  1,2, . . . , j/2-1  , 

qj  “  ’ 

f 

i  (4Jq)*1(l-n'1), 

j  -  J/2  , 

and 

p»  ■ 1  -  ji  pj  • 

J/2 

-  S  q  . 
j-l  3 

(61a) 

(61b) 


(62a) 


(62b) 


Pn  "  i  *  Z  P4  ’  s  <!.  •  (62c) 

j-l  3  j-l  3 

In  (61)  and  (62)  we  assume  J  is  an  even  integer. 

An  upper  bound  can  be  obtained  as  follows.  The  conditional  probability 
of  error  given  by  (58)  is  not  convex  in  (1  <  k  <  K,  k  ^  i).  However 
if  we  upper  bound  the  sum  of  squares  of  (53)  for  case  (i)  or  of  (55),  (53), 
and  (56)  for  case  (ii)  by  the  square  of  the  sum,  the  upper  bound  on  P 


becomes  a  convex  function  of  the  t^'s.  We  then  obtain  a  discrete  approximation 
to  the  integral  with  respect  to  t.  as  for  P  .  The  upper  bound  P  is  the  same 
for  cases  (i)  and  (ii)  and  the  distribution  of  am(k,i)  is  given  by  (59),  where 
the  pj  are  defined  by 


r  (^r^l  +  N"1)  ,  j  -  1,2 . J-l  , 

k  (1+J_1)(4q)"1(l  +  N‘1),  j  -  J 


(63) 


A  similar  approximation  and  upper  bound  can  be  obtained  for  the  sequences 
of  independent  random  elements. 

Finally  we  note  that  in  order  that  the  approximation  and  bound  presented 
in  this  section  be  tight  and  computationally  efficient  we  need  to  assume 
that  11^  ■  for  all  k  ^  i.  If  this  is  not  the  case,  we  can  still  work  with 
=»  maxC^} ,  but  the  approximation  and  the  upper  bound  obtained  above  are 
not  expected  to  be  very  tight,  so  that  it  might  be  preferable  to  work  with  the 
bounds  suggested  in  Section  III  which  are  not  affected  by  the  different  power 
levels . 

In  Table  2  the  approximation  obtained  in  this  section  is  compared  with 
the  Improved  lower  bound,  the  approximation,  and  the  upper  bound  of  Section 
III.B  for  the  first-order  Markov  hopping  patterns  and  p.^ * p.^  for  all  k. 

The  approximation  (for  both  cases  (i)  and  (ii)  and  the  bound  Py 
are  evaluated  for  J  *  4.  It  turns  out  that  they  are  rather  insensitive  to 
increases  in  J  as  long  as  J  ^  4.  Values  for  are  given  in  Table  2(a)  for  (i) 
constant  fading  and  (ii)  independent  fading.  The  notations  and 


respectively,  are  used  for  these  two  cases.  Independent-fading  turns  out  to  be 
the  most  favorable  case  although  the  difference  is  less  than  ten  percent.  Also  notice 


Table  2.  Bit  error  probability  for  nonaelective  Rayleigh  fading 


K  ■  5,  q  « 

a/N0  (dB) 

100 ,  and 

N.  -  10 

p(ii) 

A 

?U 

PA 

pu 

6 

1.64 

1.71 

1.72 

1.72 

1.74 

1.82 

(XIO  h 

8 

1.19 

1.25 

1.26 

1.27 

1.28 

1.37 

(X10  ,) 

10 

8.41 

8.95 

9.02 

9.09 

9.21 

10.14 

(X10-;) 

12 

5.80 

6.29 

6.36 

6.44 

6.54 

7.52 

(xio'p 

15 

3.28 

3.73 

3.81 

3.90 

3.97 

5.01 

(X10  ‘) 

20 

1.37 

1.78 

1.87 

1.97 

2.02 

3.10 

(X10  ,) 

QO 

0.44 

0.83 

0.92 

1.02 

1.06 

2.16 

(X10 

b) 


c) 


K  =  10,  q 

-  1000,  and  N, 
b 

-  10 

I/Nq  (dB) 

?L 

?u 

*A 

pu 

6 

8 

1.67 

1.20 

1.68 

1.22 

1.69 

1.22 

1.70 

1.24 

(xio'h 
(X10  t) 

10 

8.35 

8.51 

8.54 

8.74 

(X10  ,) 

12 

5.65 

5.79 

5.82 

6.04 

(X10  ,) 

15 

3.05 

3.18 

3.20 

3.44 

(X10,) 

20 

1.07 

1.20 

1.22 

1.46 

(xio“;> 

00 

0.10 

0.23 

0.24 

0.49 

(X10~^) 

K  *  15,  q  * 

S/NQ  (dB) 

1000,  and  N, 

D 

PL 

=  10 

pu 

PA 

P 

_ u_ 

6 

1.66 

1.69 

1.70 

1.72 

(X10  h 

8 

1.20 

1.23 

1.23 

1.26 

(xio'h 

10 

8.36 

8.60 

8.65 

8.97 

(Xio  ,) 

12 

5.67 

5.90 

5.94 

6.28 

(Xio  ,) 

15 

3.08 

3.30 

3.33 

3.69 

(xio-;) 

20 

1.12 

1.33 

1.35 

1.73 

(X10"“) 

00 

0.16 

0.36 

0.38 

0.76 

(X10'Z) 
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that  the  bound  P  (common  for  both  cases)  differs  from  P  of  (i)  or  (ii)  by 

U  A 

at  most  twenty  percent;  therefore,  in  Tables  2(b)  and  2(c)  we  present  data  on 
Py  only  (not  on  PA) •  The  purpose  is  comparison  with  PL>  PA>  and  Py.  In 
comparing  P  and  P  we  note  that  P.  appears  to  be  an  upper  bound  for  the 
nonselective  Rayleigh  case.  Also  the  results  of  Table  2(c)  show  that  for 
q  *  1000,  K  »  15,  ■  10,  and  5/Nq  ^  20  dB  the  results  of  the  Table  show  that 

Py  ^  1.17  P  ,  ?A  —  1*3  Py  and  Py  <  1.52  Py.  Similar  observations  can  be  made 
for  the  data  provided  in  Tables  2(a)  and  2(b).  As  a  final  comment  we  point  out 
that  since  the  approximations  P  and  the  bound  P  are  expected  to  be  very  close 

A  U 

to  the  true  probability  of  error,  their  favorable  comparison  with  the  simpler 

bounds  PT  and  P  and  the  approximation  P  strongly  suggests  the  use  of  the 
1j  u  a 

latter  for  the  design  of  SFH/SSMA  systems. 
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V.  SELECTIVE  FADING 

In  this  section  we  consider  a  general  wide-sense  stationary  uncorrelated- 
scattering  (WSSUS)  fading  channel.  This  model  is  described  in  detail  in  [1] 
and  [8,  Ch.  9]  and  is  employed  in  the  analysis  of  direct-sequence  SSMA 
communications  over  fading  channels  in  [4].  We  assume  that  f »  qA'»  so 
that  narrow  band  signal  models  can  be  employed.  The  input  to  the  k-th  channel 
is  sic(t“Tjc)  where 

sfe(t)  -  Re{xk(t)exp(j  2nfct  )}  (64) 

and 

^(t)  =»  v^P  exp{j(2  tt  [bk(t)A+  fk(t)  ]t  +9k(t)  +ak(t))}  .  (65) 

The  corresponding  output  is  yk(t  "Tk)  where 

yfc(t)  ■  Y0sk(t)+Re[uk(t)exp(j2nfct)}  (66) 

and 

^(t)  -  Yk  J  hk(t»T)xk(t -T)dT  ,  (67) 

so  that  the  received  signal  for  this  channel  is  given  by  (9). 

If  Yq  *  1  then  there  is  a  (non-faded)  specular  component  present  in  the 

output  of  the  channel,  and  the  channel  is  a  kician  fading  channel  (as  in  [4]). 

2  2 
In  this  case  Yk  plays  the  same  role  as  the  parameter  Y  of  Section  III.  If 

Yq  ■  0  there  is  no  specular  component,  and  the  channel  is  a  Rayleigh  fading 

2 

channel.  In  this  case  Yk  plays  the  same  role  as  the  parameter  p>k  of  Section 


IV. 
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The  fading  process  hk(t,r)  (which  can  be  thought  of  as  the  time-varying 
impulse  response  of  a  lowpass  filter)  is  a  zero-mean  complex  Gaussian  random 
process  with  autocovariance 

Eth^t.-Oh^s.a)]  -  P  k(t  -  s,t)6  (t  -ct)  ,  (68) 

where  6(0  is  the  Dirac  delta  function  and 

m 

J*  Pk(O.T)dT  »  1  . 

-OB 

Two  special  cases  of  the  model  considered  in  [1]  and  [4]  are  the  purely 
time-selective  and  purely  frequency-selective  WSSUS  fading  channels  (see 
also  [2]  and  [3]). 

In  the  present  paper  we  consider  a  somewhat  more  general  model  which  is 
both  time  and  frequency  selective-  This  is  a  special  doubly-dispersive  model 
that  is  characterized  by 


Pk(t-s,T)  *  rk(t  -  s)gfc(T)  .  (69) 

If  g^CO  s  6(t)  the  channel  is  not  frequency  selective.  If  rk(5)  ®  1  it  is 
not  time  selective. 

As  usual  ([l]-[4]),  some  limitations  are  imposed  on  the  selectivity 
of  the  channel.  First  it  is  assumed  that 

gk(T)  ^  0  for  Jt  |  >  T,  (70) 

which  is  a  constraint  on  the  frequency  selectivity  of  the  channel  that  allows 
us  to  restrict  attention  to  the  intersymbol  interference  from  the  two  adjacent 
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data  bits.  This  assumption  can  be  relaxed,  but  the  error  probability 
computations  become  more  difficult.  The  second  assumption  is  that  two 
signals  which  are  transmitted  at  different  frequencies  have  non-overlapping 
spectra  at  the  receiver.  This  is  primarily  a  limitation  on  the  time 
selectivity  of  the  channel,  but  it  also  is  related  to  the  spacing  A. 

The  analysis  of  the  receiver  follows  that  of  Section  IV,  so  many  of  the 
details  are  omitted.  The  output  of  the  in-phase  component  of  each  of  the 
two  branches  of  the  i-th  receiver  is 


Z  =  vn(D  +1  )  +  (P/8)^T(v.F  +  2  v.  +N  .  (71) 

c,m  ‘0V  c,m  c,m  Ti  c,m  ‘k  c,m  c,m 


The  terms  D  ,  I  ,  and  N  are  as  in  Section  IV  if  we  replace  A 


(i) 

c,m'  c,m'  c,m  ‘ 

by  1  and  by  0.  The  terms  F  and  I^k’^  are  (normalized)  faded 

l  c,m  c,m  v  ' 

versions  of  the  desired  signal  and  the  multiple-access  interference  due  to 
the  k-th  signal.  These  terms  are  defined  for  the  X-th  decision  bit 
(X  -  jNb  +  p)  by 


F  -  Re(F  ) 
c,m  m 


and 


f(k»i)  ,  Re(I^k,i))  , 
c,m  v  m  '  * 


(72) 


(73) 


where 


F  ■  T-1  J*  J  ht(t,T)T)t  m(t,T)exp[  j  ^  (t,T)]dTdt  (74) 

XT  ’ 

-1  n(X+1)T,* 

Im  '  *T  I  ^  hi<c-rk’T)\,i,m(t*Tk  +  T)exPl^k,i(t*Tk  +  T)ldTdt:' 


(X  +  1)T 


and 


(75) 


In  (74)  and  (75)  the  functions  h^  ^  m  and  ^  are  given  by 


i 


J 
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\,i,m(t,T)  =  6[fk(t  *T),fi(t)]5[bk(t  -r),m] 


(76) 


and 


tkji(c>T )  *-2tt  [fc  +  bk(t -T)A  +  fk(t  -  t)]t  +0k(t  -  t)  +ak(t  -  t)  -0t(t>.  (77) 


The 


functions  T) .  .  and  y  .  are  denoted  by  "*1  and  ^  ,  respectively. 

1  j  X  jIQ  1)1  i  ,m  1 


Notice  that  F  is  nonzero  if  and  only  if  both  f,  (t-T)  *  f.  (t)  and 
m  i  i 

b,  (t-T)  =  m  for  some  t  and  t  (similarly  for  I  *  ').  This  is  a  result  of 
i.  m 

our  assumption  for  the  time -selectivity  of  the  channel  and  the  size  of  A- 

In  the  analysis  below,  the  expectations  and  probabilities  are  conditioned 
(k) 

on  b^  and  rk  for  1  <  k  <  K.  However,  the  error  probabilities  that  are 

(k) 

obtained  do  not  depend  upon  Tk  or  b^  for  k  t  i.  So  in  the  last  step  we 
only  have  to  average  over  b^\ 

A.  WSSUS  Rayleigh  Fading  Model  (yq  *  0) 


The  bounds  of  (32) -(34)  and  the  approximation  of  (35)  are  employed 

except  that  fading  must  be  accounted  for  in  Pq  and  P.  For  Yq  =  0  and  K  =  1, 

Z  is  the  sum  of  two  random  variables  ((P/8)^T  F  and  N  )  of  which 
c,m  c,m  c,nr 

the  first  is  conditionally  Gaussian  and  the  second  is  Gaussian.  Furthermore, 

it  is  not  hard  to  see  that  Z  ,  and  Z,  .  are  conditionally  independent,  and 

c ,  i  c ,  — i 

so  are  Z  .  and  Z  ..  Since  a  *  a  then  [8,  p.  587] 
s,l  s ,  — 1  c,m  s,m  1  ’  r 


2  2  2  -1 
P  n  -  ct  .  (a  .  +a  .) 

3>0  c,-l  c,l  c,-l 

is  the  conditional  probability  of  error  given  there  are  no  hits  where 


(78) 


CTc,m  ■  <”2'8>n  v«!rc,»J  +  V'16 


(79) 
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It  is  convenient  to  normalize  a  and  write  (78)  as 

c  ,m 

Pe,0  =  v-l(vl  +  v-l)"1  (80) 

_  2 

where  v  is  given  in  terms  of  <S  =  Y4  PT  by 
m  1 

\  ■  2  VartFc,J  +  •  •  <«» 


The  expression  for  Var{Fc  depends  on  the  position  of  the  data  bit  within 

the  interval  [  jT^,  (j+l)T^) .  For  the  p-th  bit  of  the  j-th  hop  (X  *  jN^  +  p)  de 
define  6m  -  5(l,m),  6^  =  so/^.m),  and  6j^  «  6(b^,m).  Let 


Ht(v) 

and  define 

*  2T-2  r 

J  0  (v  -u)ri(u)du, 

T 

(82) 

Fi“* 

f  gt(T)H.(T)dT, 

0 

(83a) 

A 

T 

ft 

F.  = 

1 

and 

f  g1(T)H,  (T -T)dT  , 

0 

o  T  T  T  -T 

(83b) 

Gi  31  ' 

r  r  gt(r)  r  r  r  (t-s)dtds  . 

0  1  Vo  1 

(83c) 

The  following  expressions  for  Var[Fc  m}  are  derived  in  the  Appendix. 
First  for  p  *  0  we  find 


V“trcJ 


a»Fi 


(84) 


For  p  ■  -  1  (84)  is  valid  provided  we  replace  5^  by  6^.  Finally,  for 

0  <  p  <  -  1,  the  expression  is 


v«Cfc>J 


%[(6’ +5")F. 
©  mi 


26«ri 


26.(6' +«")G,1 
mm  m  i 


(85) 


For  the  first-order  Markov  hopping  patters  and  the  RS  hopping  patterns  the 
quantity  6q  that  appears  in  (84)  is  identically  0.  For  the  sequences  of 
independent  random  elements  6q  is  a  random  variable  with  P{6q“  1}  "  q"^  and 
p[6q  *  °}  -  1  -  q"1  . 

Notice  that  for  0  <  p  <  -  1  (i.e.  for  the  internal  bits  of  each 

dwell  interval),  Var[F£  ^  does  not  depend  on  the  hopping  pattern. 

It  turns  out  that  the  average  probability  of  error  for  these  bits 
(0  <  p  <  N^  -  1)  is  larger  than  that  of  the  first  and  last  bits  (p » 0  and 
p=N,  -  1)  .  Thus  we  use  (85),  and  not  (84),  in  order  to  obtain  an  upper  bound 
on  Pq  which  applies  for  all  values  of  p.  As  a  consequence  of  using  (85),  we 
obtain  a  bound  on  which  does  not  depend  on  the  hopping  pattern. 

In  order  to  obtain  the  limiting  error  probability  (as  the  channel  becomes 
nonselective)  it  suffices  to  let  g^ (t )  -  6(t)  and  r^u)  =  1.  We  then  have 
F.^  »  G.^  =*  0  and  F^  -  %  so  that  PQ  is  given  by  (41)  with  A  *  5/Nq  »  y^  p  . 
Similarly,  to  obtain  the  irreducible  error  probability  (as  p  —  •)  we  simply 
disregard  the  second  term  in  the  right-hand  side  of  (81) . 

For  the  WSSUS  Rayleigh  fading  model  we  say  that  a  hit  occurs  from  the 
k-th  signal  whenever  t.  ,  b,  (t),  and  f.  (t)  are  such  that  Var[l^’^^J  #  0. 

The  probability  P  of  such  a  hit  depends  upon  N^  and  q.  For  the  first- 
order  Markov  hopping  patterns  we  have 

,J£P(j4i  .  (86, 

b 

In  deriving  (86)  we  used  the  fact  that  for  the  selective  fading  model  used 
in  this  section,  as  many  as  4  adjacent  bits  from  the  k-th  signal  may 


interfere  with  each  bit  of  the  i-th  signal.  The  expressions  (25)  and 
(32)- (35)  apply  with  P  replaced  by  P^  and  Pq  evaluated  as  explained  above. 
For  memoryless  hopping  patterns  the  corresponding  result  is 


p-pu4?ii  +  ^<1-5>i  •  <87> 

Both  bounds  in  (86)  and  (87)  are  tight  for  N^  a  3 .  For  the  RS  hopping 
patterns  the  corresponding  result  is 


p  <  P 

u 


(88) 


Notice  that  the  bound  in  (88)  is  the  same  as  in  (23)  which  was  obtained  under 
nonselective  fading  conditions.  This  is  due  to  the  fact  that  the  RS  hopping 
patterns  do  not  repeat  within  a  period. 

B.  WSSUS  Rician  Fading  Model  (yQ  »  1) 

In  this  case  the  conditional  error  probability  given  there  are  no  hits 
is  [8,  p.  587] 


Pe,0  - 


(89) 


Upon  normalization,  (89)  reduces  to 


e,0 


v-l(vl  +  v-l 


)_1  expC-[Y^(v1  +  v_l)r1} 


(90) 


where  v  is  defined  by 
m 

vm  -  2  Var{Fc>m}  +  (I  +  Y^KyJ  */Nq)"\  (91) 

<J/Nft  »  (1  +  Y.,)P»  and  Var^F  }  is  given  by  (82)-(85).  Finally  in  order  to 
u  i  c,m 

obtain  Pq  we  have  to  average  Pg  q  with  respect  to  the  data  bits  (^x+l’^X-l^ * 
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For  Rician  fading  the  hits  from  the  k-th  signal  may  occur  from  either 
the  direct-path  component  or  the  faded  component.  The  probability  of  a  hit 
from  the  k-th  signal  due  to  the  direct-path  component  is  the  same  as  for 
nonselective  Rayleigh  fading  (this  was  evaluated  in  Section  III).  The 
probability  of  a  hit  due  to  the  faded  component  is  evaluated  above  (for 
Rayleigh  fading) .  The  union  bound  provides  simple  upper  bound  on  the 
probability  of  a  hit.  This  is  given  by 

p  £  p;  - 1  d +  h  <92) 

b 

for  first-order  Markov  hopping  patterns  and 

p  f1  +  ^-n)}  (93) 

u  q  Nb  q 

A 

for  memoryless  random  hopping  patterns.  For  RS  hopping  patterns  P  is 
still  bounded  as  in  (88);  that  is, 

P  <P^~PU  •  (94) 

By  substituting  for  Pq  in  (32)-(35)  and  replacing  P  by  P^  in  (25)  we  have 
lower  bounds,  an  approximation  and  an  upper  bound  on  the  average  probability 
of  error. 

In  Tables  3  and  4  the  approximation  given  in  (35)  is  obtained  for 
purely  frequency-selective  Rayleigh  and  Rician  fading  channels,  respectively. 
The  system  parameters  are  K  ■  15,  q  ■  1000,  and  ■  10.  First-order  Markov 
hopping  patterns  are  employed.  The  covariance  function  of  the  frequency- 
selective  channel  is  triangular,  so  that  the  rms  multimath  spread  o  defined 

2 

by  a  ■  J  t  g(T)dT  is  related  to  the  parameter  d  of  [3]  by  d  ■  2.22  o/T. 

mm 

We  let  “  Y  for  all  k.  Then  in  Table  3,PA  is  given  as  a  function  of 


Bit  error 

probability 

for  Rayleigh 

frequency- 

■selective 

fading 

(K  -  15, 

q  *  1000,  and 

"t  ■  10>- 

£/Nq  (dB) 

O-0.05T 

fl-0.1T 

O-0.15T 

O-0.2T 

6 

1.75 

1.82 

1.91 

2.01 

(xio'b 

8 

1.28 

1.35 

1.44 

1.54 

(xio  :) 

10 

0.91 

0.97 

1.06 

1.17 

(xiop 

12 

6.31 

6.88 

7.71 

8.84 

(xio  ,) 

15 

3.63 

4.13 

4.94 

6.08 

(xio'p 

20 

1.59 

2.04 

2.82 

3.95 

(Xio  ,) 

OO 

0.58 

1.00 

1.76 

2.89 

(X10‘Z) 

Table  4.  Bit  error  probability  for  Rician  frequency-selective  fading 
(K  *  15,  q  «  1000,  Nb  -  10,  and  o  *  0.05T). 
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Y2-l 
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8 
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10 
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12 
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15 
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(xio 

20 
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—  2  2 

5/Nq  *  (1  +  Y  )P  f°r  a  *  0.05  T  and  for  five  different  values  of  Y  •  Notice 

2 

that  as  Y  "*  •  the  probability  is  not  the  same  as  the  second  column  of 

Table  3.  Although  P^  is  the  same  in  this  limiting  case,  the  fact  that 
P^  <  P^  (compare  (86)  to  (92))  implies  that  the  two  cases  give  different 
values  of  the  bit  error  probability. 

Finally  we  compare  P^  for  nonselective  and  frequency-selective  Rayleigh 
fading  for  K  »  15,  q  *  1000,  and  ■  10,  (first-order  Markov  hopping  patterns 
are  employed).  From  Tables  2(c)  and  3  we  see  that  the  probability  of  error 
for  the  frequency-selective  case  is,  for  5/Nq  ■  12  dB  and  a  *  0.05,  1.1  times 
that  for  nonselective  fading,  and  it  becomes  1.5  times  the  corresponding 
probability  for  nonselective  fading  as  a  increases  to  0.2  T.  Similarly 
for  <?/Nq  *  20  dB  the  ratio  of  the  two  probabilities  ranges  from  1.2  for 
a  »  0.05  T  to  2.9  for  a  »  0.2  T. 


APPENDIX 


In  chi s  appendix  we  develop  Che  expressions  for  Var[Fc  ffl) .  As  in  [4] 

we  can  write  VarCF  }  as 
c  ,m 

VariF^  J  =■  E[ReCFj]2  =■  h  E[F  F*]  ,  (A-l) 

c,m  qi  m  m 

where  we  used  the  fact  that  [1]  E{h^(t,T )h^(s ,<?)}»  0.  Upon  substitution  for 
(74),  (68)  and  (69)  in  (A-l)  we  find 

2  -  (X+l)  (X+1)T 

v.rtFc>nl-i  ;XI  't<'->V,<e.T>nlja<..T) 

exp[j[*1(t,T)  -tt(s,T)]}dtdsdT  .  (A-2) 

Notice  that  T).  (t,T)Tl,  (s,r)  4  0  only  for  those  t,  s,  and  t  for  which  the 
i  ,m  l  ,m 

following  three  conditions  hold:  f^t-T)  -  ft(t),  f^s -t)  «  f^s),  and 
b^(t  - T )  »  (s  - t )  *  m.  But  these  three  conditions  imply  ( t -t)  *  at(t), 
Qfi(s  -  t)  »cr^(s),  ©^(t-T)  *9^(s-t),  respectively.  Also  ar^(t)  *  a^(s)  » 
and  0^(t)  *  P±(s)  *  f3j^  for  t  and  s  in  [XT,(X+1)T).  Consequently, 
ti(t,T)  *  tjr ^ (s ,T )  for  these  values  of  t,  s,  and  t.  As  a  result  we  may  let 

exp(j[*1(t,T)  -  tt(s,T)]3  -  1 

in  equation  (A-2) . 

The  next  step  is  to  write  (A-2)  as 

®  p-1  Nb"1 

Var(F  )  -  Z  [d(i)  Z  A  (A,n)  +  A  (i,p)  +  d(i+l)  L  A(A,n)l, 

C’m  X—  n-0  ”  “  n-|H-l  m 


(A-3) 


A- 2 


where  for  n  i  p 


A  (*,n) 
m 


fT 

gj (t+£T. +nT)[  A  (Z,a+l)F(r)  +  A  (l, n)F(r)  + 

i  n  m  in 

0 

2A  U,n+1)A  (£,n)G(T)]  dx , 
m  m 


(A-4a) 


and  for  n  ■  p 


A  (£,p) 
m 


8i(T+£VpT)[d(i+1)V*’P+1)F(T)  +  d(^)Am(£,P)F(T)  + 


2d(i+l)d(£)A  ( i , p+1 ) A  U,p)G(t)]dT. 
m  m 


(A-4b) 


In  (A-3)  -  (A-4)  we  also  need  the  definitions 


dU)  -  «(f£>.  fj(i)),  (A-5) 

V*'">  *  S<bA(\-„-  »>•  <A'b> 

b 


and  (cf.  (82)  -  (83)) 


_2 

(•XT+r 

(•XT+t 

F(t) 

*  T  ^ 

r  (t-s)dtds  -  H. (t). 

(A-7a) 

* 

XT  - 

1  XT  1  1 

r(X+l)T  f(X+l)T 

F(t) 

-  T  2 

r.(t-s)dtds  -  H. (T-t) 

(A-7b) 

* 

'XT+t 

>  XT+t  1  1 

-2 

rXT+T 

r(X+l)T  fT  rT 

G(r) 

-  T  2 

r.(t-s)dtds  *  r.(t-s)dtds. 

(A-7c) 

XT 

XT+t  1  >o  1 

Notice  that  the  result  of  (A-3)  is  quite  general  and  it  accounts  for  the 
Intersymbol  interference  due  to  many  data  bits.  However,  because  of  the 
assumption  (70)  only  the  terms  £• 0,  n-0  and  £— 1,  n-N^-1  of  (A-3)  give 
nonzero  contributions,  and  thus  (A-3)  reduces  to  (84)  -  (85). 
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APPENDIX  B 


ANALYSIS  OF  A  SLOW  FREQUENCY-HOPPED  SYSTEM  WITH  POISSON  TRAFFIC 

In  this  section  the  packet  error  probability  and  the  throughput  rate 
are  determined  for  a  particular  frequency-hopped  system  when  the  number 
of  packet  transmissions  in  a  slot  is  given  by  a  Poisson  random  variable. 
The  packet  error  rate  and  throughput  under  the  Poisson  traffic  assumption 
are  significant  in  view  of  the  local  Poisson  approximation  and  recursive 
retransmission  control  strategies  discussed  in  the  next  section. 

The  system  of  interest  is  assumed  to  be  packet-synchronized.  A 
sufficient  time-guard-band  must  be  maintained  between  packet  slots  to 
maintain  sychronization  in  the  face  of  differential  delays  due  to  the 
spatial  distribution  of  the  stations.  Synchronization  at  the  level  of 
bits  or  bytes  is  not  assumed. 

Each  packet  transmission  is  declared  successful  or  not  according  to 
some  criteria  (a  specific  choice  is  given  below) .  The  following  defini¬ 
tions  will  be  used 

r(m|k)  »  P[m  successful|k  packets  transmitted  in  slot] 
k 

r(k)  *  2  m  r(m|k) 
m-0 

and 

P(k)  -  1  -  r(k)/k. 

Thus  r(*|k)  is  the  distribution  of  the  number  of  successes,  r(k)  is  the 
mean  number  of  successes,  and  P(k)  is  the  average  probability  of  failure 
for  a  typical  packet,  all  given  that  k  packets  are  transmitted  in  the  slot. 

Similarly,  define 

rp(m|G)  -  E[r(m| K)  ] 


and 


rp(G)  =  E[r(K) ] 


Pp(G)  -  l-rp(G)/G  j 

’J 

j 

\ 

where  K  is  a  Poisson  random  variable  with  mean  G.  Thus  rp ( • | G)  and  rp(G) 
are  the  distribution  and  mean  of  the  number  of  successes  and  Pp(G)  is  the 
probability  of  failure  of  a  typical  packet,  all  given  that  the  number  of 
transmissions  in  the  slot  is  a  Poisson  random  variable  with  mean  G. 

The  specific  FH  system  will  now  be  described.  The  frequency  spectrum 
is  divided  into  q  frequency  slots  and  the  packets  are  divided  into  n  bytes  each 
Each  byte  is  transmitted  at  a  frequency  chosen  from  the  q  frequencies  with 
equal  probability,  independently  of  the  frequencies  chosen  for  other  bytes. 

It  is  then  appropriate  to  use  a  burst-error  correcting  code  —  we  will 
assume  that  a  Reed-Solomon  code  is  used.  We  will  also  assume  that  a 
packet  consists  of  exactly  one  codeword  from  a  RS  code  for  which  up  to  t 
byte  errors  can  be  corrected.  This  provides  us  with  a  natural  definition 
of  a  successfully  transmitted  packet.  A  packet  is  declared  successfully 
transmitted  if  at  most  t  byte  errors  occur.  Both  the  (31,15)-code  (with 
n»31  bytes,  five  bits  per  byte,  15  information  bytes  and  t  ■  8)  and  the 
(255,127)-code  (with  n  ■  255  bytes,  eight  bits  per  byte,  127  information 
bytes  and  t  •  64)  will  be  considered. 

In  the  following  let  X»G/q,  so  that  X  is  the  traffic  intensity  per 
frequency  slot.  Also,  P(k,q)  and  Pp(X,q)  will  be  written  in  place  of 
P(k)  and  Pp(G)  in  order  to  make  the  dependence  on  q  explicit.  Finally, 
let  q(X,q)  ■  X  (l-Pp(X  ,q)) .  Thus  n  denotes  the  average  throughput  per 
frequency  slot. 


Assume  now  that  byte  errors  are  independent  in  the  absence  of  multi¬ 
access  interference.  This  independence  assumption  is  true  for  an  additive 
white  Gaussian  noise  (AWGN)  channel.  The  assumption  is  also  approximately 
true  for  an  AWGN  channel  with  fading  if  q  is  so  large  that  very  few  fre¬ 
quency  slots  are  hit  by  more  than  one  byte  for  any  packet,  or  if  the  fading 
process  of  the  channel  model  has  a  short  correlation  time  compared  to  the 
typical  elapsed  time  between  visits  to  a  given  frequency.  Also,  by  inter¬ 
leaving  codewords,  it  is  possible  to  approximately  achieve  the  situation 
with  independent  byte  errors  even  for  relatively  slowly  fading  channels. 

Now  let  p^  be  the  byte  error  probability  in  the  absence  of  multi¬ 
access  interference.  Then  the  byte  error  probability  given  that  k 
packets  are  transmitted  in  a  slot  is 


■k 


1  -  (1  -  (f  - 

q  q 


(k  >  1) 


If  byte  synchronization  were  possible,  the  term  - 2  this  expression 

q  q 

could  be  replaced  by  1/q.  By  the  assumed  independence  of  byte  errors  in 
the  absence  of  multi-access  interference  and  the  memorilessness  of  the 
random  hopping  pattern,  the  byte  errors  (including  multi-access  interference) 
are  conditionally  independent  given  k.  Thus  the  packet  error  probability  is 


P(k,q)  -  2  (J)  p*  (1-p.) 

i-t+1 


n-i 


which  can  be  used  to  compute 


Pp(X,q)  «  £ 

k-1 


-Xq  _  s  k 

* - k  P(fc,q)Aq 


and 


n(X,q)  -  X(l-Pp(X,q)) 

which  are  the  desired  packet  error  probability  and  throughput  per  frequency 
slot  for  Poisson  traffic. 

The  computation  of  Pp(X,q)  (and  hence  also  Tj(X,q))  simplifies  in  two 
special  cases:  First,  when  q*l, 

Pp (X ,  1)  -  l-e_X  (1  -  P(l,l)) 

since  P(k,l)  *1  for  k  >  2.  This  implies  that 

•7 (X ,  1)  -  *e"X  (1-P(1,1>). 

That  is,  when  q*l,  the  throughput  is  Xe  X  (which  is  the  throughput  for  a 
noiseless  slotted-ALOHA  channel)  times  the  success  rate  in  the  absence  of 
multiple-access  interference. 

The  second  special  case  is  obtained  by  letting  q  and  G  tend  to  infinity 
with  X  * G/q  fixed.  The  limiting  packet  error  probability  is  then 


P  (X,+  ~)  -  lim  Pp(X,q)  -  lim  P(k,q) 

*■  q  -*■  oo  *■  q  ->•«> 

k  -*-oc 
X-k/q 

-  *  s  (?)  pi  (i  -  pj11'1 

i-t+1 

where 

Poo  •  li®  Pk  ■  1  -  e~2X  d  -  Pi)  • 
q  -*-»o 

k 

X-k/q 


Numerical  results  are  given  in  Figs.  1-4  and  in  Tables  1-3.  We  see 
in  Fig,  1  that  for  no  channel  noise  and  using  the  (255,127)  Reed-Solomon 
code,  a  smaller  packet  error  probability  is  achieved  by  q  =  +°°  than  by 
q *  1  if  and  only  if  X  is  smaller  than  about  0.13.  For  an  intuitive 
understanding  of  this  it  is  important  to  keep  in  mind  the  following  two 
facts.  First,  at  the  level  of  byte  errors,  the  essential  effect  of  varying 
the  parameter  q  is  that  as  q  decreases,  the  occurrences  of  byte  errors 
within  a  single  packet  become  more  positively  correlated.  Secondly,  since 
X  is  the  traffic  normalized  per  frequency  slot,  for  fixed  X  the  byte  error 
probability  and  therefore  also  the  mean  number  of  byte  errors  per  packet 
does  not  strongly  depend  on  q.  Summarizing  these  two  facts,  for  larger  q 
the  distribution  of  the  number  of  byte  errors  tends  to  be  more  tightly 
concentrated  near  the  (almost  q-independent)  mean  number  of  byte  errors. 

Thus,  whether  or  not  the  packet  error  probability  is  smaller  for  large 
q  than  for  small  q  is  determined  by  whether  or  not  the  error  correcting 
capability  of  the  code  can  accomodate  any  number  of  byte  errors  "near" 
the  mean  number  of  byte  errors.  Since  the  mean  number  of  byte  errors 
increases  with  X ,  it  follows  that  for  small  enough  X  the  packet  error 
probability  is  smaller  for  large  q,  and  conversely  for  large  X  the  packet 
error  probability  is  smaller  for  small  q. 

When  the  byte  error  probability  in  the  absence  of  multi-access 
interference  p^  is  increased  from  zero  to  0.1,  the  packet  error  probability 
does  not  significantly  increase  for  q  ■  1  while  it  does  for  larger  values 
of  q.  (Compare  Figs.  1  and  2.)  As  a  result,  the  value  of  X  at  which  the 
packet  error  probability  for  q-°°  surpasses  the  packet  error  probability 
decreases  to  X  -0.078.  (See  Fig.  2.)  Thus  in  the  presence  of  channel  noise, 
the  crossover  value  of  X  can  become  quite  small. 


Turning  to  Figs.  3  and  4  we  observe  that  the  maximum  throughput  (over 
all  X)  is  much  greater  for  q  =  l  than  for  q*°°.  However,  the  maximum 
throughput  for  q  = 1  can  only  be  achieved  by  maintaining  a  mean  traffic 
intensity  X*1  which  causes  the  packet  error  probability  to  exceed  .63. 
Hence,  although  q=l  offers  greatly  increased  maximum  throughput,  the 
increase  comes  at  the  expense  of  either  many  retransmissions  (which,  if 
possible  at  all,  generally  increase  delays)  or  a  high  packet  loss  rate. 

Discussion  of  Method 

The  method  of  using  the  local  Poisson  approximation  as  discussed  here 
and  in  the  next  section  is  admittedly  only  an  approximation.  It  is  impor¬ 
tant  to  emphasize  however  that,  as  shown  in  [7],  the  method  does  lead  to 
channel  stability  (even  taking  approximations  into  account) . 

Another  approach  to  the  analysis  of  delay  in  a  random-access  system 
would  be  to  use  a  more  detailed  model  of  the  transmitters  —  allowing  them 
to  obtain  multiple  packets  and  then  buffer  delay  could  be  discussed.  For 
such  analysis  so  far  in  the  literature,  the  total  system  is  usually 
described  as  a  (many  state)  Markov  chain.  For  such  analysis,  the  main 
obstacle  has  been  the  large  size  of  the  state-space.  Here  we  wish  to 
point  out  another  difficulty  which  arises  for  such  detailed  analysis  when 
one  considers  spread-spectrum  systems.  The  problem  is  that  a  detailed 
exact  analysis  would  require  knowledge  of  the  conditional  distribution 
r('|k)  of  the  number  of  successes  given  k  transmissions  (whereas  our 
analysis  only  required  use  of  the  mean  number  of  successes) .  Some  authors 
propose  (implicitly)  that  the  distribution  of  the  number  of  successes  is 
binomial  under  the  assumption  that  the  outcomes  of  transmissions  of  distinct 
packets  form  independent  events.  It  is  clear,  due  to  the  mutually  destruc¬ 
tive  effect  of  collisions  that  this  assumption  is  not  true. 


In  summary  —  before  more  detailed  models  can  be  effectively  used, 
the  distribution  r(*|k)  must  be  better  characterized.  Nevertheless,  we 
have  found  retransmission  control  schemes  which  insure  stable  throughput 
even  without  knowledge  of  this  distribution  (see  next  section). 


Table  1.  Packet  error  probability  and  throughput  vs.  traf f ic-intensity-per 
frequency-slot  X  for  q  =  1  (no  hopping  during  packet  transmission) 
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Table  3.  Byte  error  probability,  packet  error  probability  and  throughput  vs.  traffic-intensity-per-f requency- 
slot  for  q  =  °°. 
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Throughput  vs.  traffic-intensity-per-frequency-slot  X.  Independent 
byte  errors  —  error  probability  p^“0.1. 


APPENDIX  C 


RECURSIVE  RETRANSMISSION  CONTROL  —  APPLICATION  TO 
FREQUENCY-HOPPED  (FH)  SPREAD-SPECTRUM  SYSTEMS 

Although  the  models  of  the  user  population  and  the  feedback  Information 
are  quite  simple,  the  concepts  of  f 5 J  — f 7 ]  readily  extend  to  more  complex  and 
realistic  settings.  In  order  to  illustrate  this  point,  we  shall  briefly 
describe  how  the  decentralized  dynamic  control  procedure  in  [5]  can  be  adapted 
to  the  frequency-hopping  system  described  in  the  previous  appendix.  For 
definiteness,  suppose  that  the  Poisson  model  in  [5]  is  used  to  describe  how 
the  population  of  stations  acquires  packets  to  be  transmitted. 

Our  research  has  shown  that  it  is  desirable  for  the  users  to  have  some 
feedback  information  in  order  to  suitably  control  the  traffic  level.  We 
shall  now  describe  a  method  for  the  users  to  obtain  such  Information  which 
is  appropriate  for  use  in  the  Navy's  ITF  network.  During  each  slot,  user  a 
uses  a  random  hopping  pattern  to  hop  a  receiver  among  the  q  frequencies. 

The  pattern  has  the  same  distribution  as  patterns  used  to  transmit  packets. 

For  each  dwell  time  the  user  decides  (by  a  simple  threshold  test)  whether 
or  not  the  channel  was  free  during  that  dwell  time  in  the  frequency  monitored. 
The  user  then  simply  counts  the  number  of  dwell  times  in  the  slot  for  which 
it  was  decided  that  the  channel  frequency  was  not  free.  Let  Y®  denote  the 
count  of  user  a  for  slot  t.  The  variables  Y^  comprise  the  feedback  infor¬ 
mation  upon  which  the  retransmission  control  strategy  described  next  is  based. 

Following  [5],  we  suggest  that  user  a  recursively  computes  the  sequence 
f®  via  the  multiplicative  rule 

f®+1  -  min(f“a(Yj),  1)  (1) 

Then  if  user  a  has  a  packet  to  transmit,  it  transmits  it  in  slot  t  with 
probability  f®.  A  possible  choice  for  the  function  a  is 


(2) 


a(y)  -  1  -  -  p*)r 

where 

■k  k 

p  =  1  -  exp (-2X  ) 
k 

and  X  is  a  desired  value  of  the  traffic  intensity  per  frequency  slot. 

To  understand  this  choice  of  retransmission  policy,  note  first  that 
if  the  number  of  packet  transmissions  in  a  slot  is  Poisson  with  mean 
G«qX,  then  the  probability  that  a  frequency  slot  is  used  during  any  portion 
of  a  given  one-byte  dwell  time  is 


p  *  1  -  exp (-2X ) 

k  k  k 

When  X«X  ,  p  -  p  .  That  is,  p  is  the  probability  that  a  frequency  is 

occupied  during  a  given  one-byte  dwell  time  when  the  number  of  packet 

* 

transmissions  is  Poisson  with  the  desired  mean-per-frequency-slot  X  . 

Now  given  the  current  values  of  the  retransmission  probabilities  f^ 
for  all  a  and  given  the  set  Afc  of  users  which  have  a  packet  to  transmit, 
the  conditional  distribution  of  the  number  of  packet  transmissions  in  slot 
t  is  approximately  Poisson  with  mean 


2 

a€ 


f“  +  H 


(4) 


(where  p  is  the  rate  at  which  new  packets  are  transmitted)  by  the  local 
Poisson  approximation  described  in  [5]  and  [9].  Hence, 


E  [Y“|f“,  all  a]  -  qpt 

where 

Pt  -  1  -  exp(-2Gt/q) 


(5) 


Thus  by  (1)  and  (4) 


E  [Ct+ilf“.  aii  a]  =  E  [  2  f“a(Y“)  +  /t|f“,  all  a] 


a  6  A 


t+1 


»+  2  f®E[a(Y®)|f“  all  o] 

*€Vl 


and  since  (by  (2)  and  (5)) 

Ya  A 

E  [a(Y®)If“,  all  a]  *  E  [1-7  (-jp  -  p  }  |f®,  all  a] 

=  1  -  Y(Pt  -P*) 

*  l+7(exp(-2Gt/q)  -  exp(-2G*/q)) 

we  have 

E  a11  °1  a  {l+7(exp(-2Gt/q)  -  exp(-2G*/q))}Gt 


Hence,  Gt+1  will  tend  to  be  larger  than  Gt  (resp.  smaller  than  Gt>  if  Gt  is 

* 

smaller  (resp.  larger)  than  the  desired  traffic  level  G  .  That  is,  Gfc  will 

* 

tend  to  drift  toward  the  desired  traffic  level  G  . 

Using  the  methods  of  [  7]  it  can  be  shown  that  with  this  transmission 

policy  the  channel  is  stable  whenever  the  input  rate  ju  is  smaller  than  the 

* 

average  throughput  rate  qv (X  ,q)  corresponding  to  the  desired  traffic 

* 

intensity  per  frequency  slot  X  . 


APPENDIX  D 


ACKNOWLEDGEMENT  BASED  RETRANSMISSION  CONTROL 

In  [9]  we  define  and  analyze  the  class  of  Acknowledgement  Based  Retrans¬ 
mission  Control  (ABRC)  schemes  for  random  access.  Such  schemes  require  no  more 
feedback  information  than  the  original  ALOHA  scheme  in  that  each  user  need 
only  learn  whether  or  not  its  own  transmission  is  successful.  Such  small 
feedback  requirements  are  desirable  in  a  spread  spectrum  environment  when 
channel  monitoring  is  difficult.  On  the  basis  of  approximations  developed 
in  [5]  under  this  contract  and  an  equilibrium  analysis,  we  have  found  that 
such  schemes  can  provide  satisfactory  performance  for  both  an  infinite  and 
finite  number  of  users,  as  long  as  the  retransmission  probabilities  are 
properly  chosen. 

We  also  introduced  the  possibility  that  after  a  packet  has  collided  a 
certain  number  of  times  then  it  is  rejected  and  no  longer  retransmitted. 

It  appears  that  the  probability  that  a  given  packet  must  be  rejected  can 
be  kept  to  a  satisfactorily  low  level,  while  allowing  rejections  improves 
stability  considerably.  In  Fig.  1  the  average  probability  of  success  and 
the  probability  of  ultimate  rejection  P  are  given  as  a  function  of  the 
allowed  number  of  retransmissions  k  for  the  infinite  user  model  with  input 
rate  \»0.3.  For  k  *  12  the  probability  of  rejection  is  about  one  in  ten 
thousand.  However  if  k  is  chosen  to  be  too  large  then  undesirable  bistable 
behavior  appears. 
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APPENDIX  E 


THE  METHOD  OF  MARKOV  PROCESSES  WITH  PHASES 

Another  area  of  progress  under  this  contract  has  been  the  development 
of  a  numerical  technique,  the  "method  of  phases",  which  we  have  discovered 
is  suitable  for  evaluating  certain  random  access  algorithms  in  the  presence 
of  fluctuating  traffic  rates.  Our  main  motivation  is  that  the  usual 
Bernoulli  or  Poisson  models  of  arrival  processes  are  not  "bursty"  enough 
to  realistically  model  traffic  which  random  access  schemes  are  likely  to 
face  in  practice.  So  far  the  method  has  been  successful  for  evaluation  of 
TDMA  with  buffered  users  and  varying  arrival  rates.  The  method  will  now 
be  briefly  described,  following  [8]. 

Consider  a  TDMA  system  with  m  users  and  arrival  rate  a  packets/slot/ 
user.  Let  Nt  denote  the  number  of  packets  in  the  first  user's  buffer  and 
let  9t  €  {l,...,m}  denote  which  user  is  transmitting  during  slot  t.  (9  is 
the  "phase"  of  the  system  relative  to  the  first  user.)  Then  (Nt,9t>  can 
be  modelled  as  a  discrete-time  Markov  chain  on  Z+  x  [l,...,m]  with  transi¬ 
tion  matrix  (in  the  following  matrices,  only  non-zero  entries  are 
indicated) : 

A00  A01 

A10  A1  A2 

Aq  Al  A2 

A0  A1  *2 

•  •  • 

•  •  • 

•  • 

where  the  blocks  of  P  are  the  mXm  matrices 
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In  [8]  a  general  method  Is  presented  for  obtaining  the  invariant  distribu¬ 
tion  of  such  chains.  The  general  procedure  is  as  follows  [8]: 

Step  1:  Compute  B,  where  B  is  the  minimum  non-negative  solution  to 
the  equation 


B  =  BA2BA0(I-A1)'1+  (I-Aj)"1 

Successive  substitutions  starting  with  B ^  =  0  yields  an  increasing 
sequence  converging  to  B. 


Step  2;  ^Find  tTq,  the  invariant  distribution  for  the  mXm  transition 
matrix  P*  *  AQ0°+  A^HAy,. 

Step  3:  Compute  the  constant 


c  -  TT0(I+A01B(I-R)_1)e 


where  e  *  (1,1 . 1;  and  R  ■  AjB. 

•fa 

Then  the  invariant  distribution  for  P  is  x  *  (x^,x^,...)  where 


and 


Vc 


xoaoibr' 


k-1 


k  a  1 


In  [8]  analogous  results  are  also  derived  when  P  is  truncated  to  a  finite 
number  of  levels  and  boundary  states  are  added.  This  provides  a  computa¬ 
tionally  tractable  method  to  analyze  queues  with  finite  buffers. 

Using  this  approach,  the  invariant  distribution  for  the  TDMA  example  has 
been  found  analytically  and  numerically  [10].  Our  results,  such  as  expres¬ 
sions  for  the  average  backlog,  agree  with  those  obtained  by  other  methods. 

The  advantage  of  this  approach  is  that  it  readily  extends  to  the  case  when 
the  arrival  rate  fluctuates  according  to  an  tinder lying  Markov  process,  for 
then  the  system  still  has  the  same  form  as  above  but  for  different  choices 
of  the  A^'s.  This  extension  has  been  carried  out  and  is  presented  in  [10]. 

An  example  of  our  numerical  results  are  presented  in  Fig  2.  In 
this  example  the  number  of  users  M  was  taken  to  be  4  or  10.  In  each  case 
the  mean  arrival  rate  was  p  =»  96%  of  the  channel  capacity,  and  the  actual 
arrival  rate  fluctuated  between  two  different  arrival  rates,  where  the 
switching  was  governed  by  a  two  state  Markov  chain.  The  dashed  lines 
correspond  to  an  example  when  both  of  the  input  rates  were  chosen  below 
the  system  capacity  while  the  solid  lines  correspond  to  an  example  when  one 
of  the  two  rates  is  above  the  system  capacity.  The  curves  give  the  average 
backlog  N  for  a  given  user  versus  y,  where  y  is  a  parameter  which  indicates 
how  fast  the  rate  is  switched.  For  small  y  the  switching  processes  is 
slow  so  that  a  large  backlog  results  when  one  of  the  two  rates  is  above  the 
system  capacity  (see  solid  lines). 

Our  reason  for  studying  arrival  processes  with  varying  rates  is  that 
we  feel  it  provides  a  more  realistic  model  of  bursty  traffic  than  does  the 
usual  Poisson  arrival  model.  We  are  now  in  the  process  of  analyzing  other 
random  access  disciplines  in  the  presence  of  varying  traffic  rates.  We 
suspect  that  many  random  access  schemes  will  perform  more  favorably  relative 
to  the  performance  of  TDMA  when  the  traffic  arrival  rates  vary  dynamically. 


